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I. INTRODUCTION 



Conformal Field Theory (CFT) is a powerful approach which determines the low-energy prop- 
erties of one-dimensional (ID) quantum systems or two-dimensional classical systems close to criti- 
cality [1, 2]. In particular, it allows for a complete classification of ID gapless phases and quantum 
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critical points with Lorentz invariance [3]. The low-energy gapless spectrum of the corresponding 
ID lattice model with continuous symmetry, like the spin-1/2 Heisenberg chain, is described in 
terms of representations of certain Kac-Moody (KM) current algebra [4]. This affine symmetry 
fixes the operator content of the system and the scaling dimensions of the operators are in turn 
determined by the conformal invariance of the underlying Wess-Zumino-Novikov-Witten (WZNW) 
model built from the currents of the affine symmetry [2]. The ID gapless phases are then labelled 
by the central charge c of the corresponding WZNW CFT. The simplest example of this description 
is the Luttinger liquid universality class which accounts for the low-energy properties of ID metals 
[5-7]. Such a CFT has a U(l) affine KM symmetry with central charge c = 1. The critical theory 
of half-integer Heisenberg spin chains corresponds to the SU(2)i WZNW whereas the family of 
integrable spin chains of arbitrary spins 5 admits for a critical theory the SU(2)2s WZNW model 
with central charge c = 35/(5 + 1) [8]. 

In sharp contrast to gapless systems, a general classification of all possible zero temperature 
ID gapped phases is still lacking despite numerous works in the past two decades. To the best of 
our knowledge, there exists no general symmetry based arguments that allow for a quantitative 
description of ID quantum problems with a spectral gap. A natural starting approach is to identify 
the massive phases that occur in the vicinity of quantum critical points or close to gapless phases 
which are described by some known CFTs. Very special perturbations could be singled out which 
correspond to massive integrable deformations of the original CFT [9]. The integr ability provides, 
in turn, a non-perturbative basis to construct the massive quasiparticles defining the low-energy 
spectrum of the fully gapped phase. The stability of these excitations can then be investigated 
perturbatively with respect to other non-integrable generic perturbations [10] and it will define the 
extension of the massive phase. A second more phenomenological approach relies directly on orders 
by considering very special ones and the use of extended symmetries which unify them. A typical 
example is the SO (5) theory for the competition between cl-wave superconductivity and antiferro- 
magnetism, where the order parameters of <i-wave superconductivity and antiferromagnetism are 
combined to form a unified order parameter quintet [11]. Several different competing orders may 
be transformed into each other under this extended global symmetry. 

In this work, we will combine these two approaches to describe spin-gapped phases for a large 
class of ID fermionic models where the gap is exponentially small in the weak coupling regime. 
More precisely, we will be concerned with weakly interacting iV-component fermionic models such 
as n-leg electronic spin-1/2 ladders (N = 2n) [6, 7] or one-dimensional fermionic cold atoms with 
hyperfine spin F (N = 2F + 1) [12, 13]. In absence of interactions, those systems display a large 
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degeneracy, with a corresponding large symmetry group mixing microscopic degrees of freedom, 
to be called G in the following . For the aforementioned models, the group G will be SU(iV) or 
SO(2iV), depending on the filling. This degeneracy allows for many possible competing orders. 
The low-energy properties of the non-interacting system are governed by a CFT built on the 
group G: the G& WZNW CFT. For the applications to coupled fermionic chains, the level of this 
CFT is k = 1, but most of our analysis can be carried out in the general k case. Interactions 
generically break the symmetry G down to a subgroup H which is the "true" physical symmetry 
of the problem. The degeneracy lifting which arises from the interactions results in general in the 
opening of a spectral gap, which is accompanied by the stabilization of one of the possible orders. 
Questions we will address in the following are: what kind of competing orders can be stabilized? 
How do they depend on the algebraic properties of the groups involved in the symmetry breaking 
scheme G H? Can these orders be related between themselves? And what is the nature of the 
(zero-temperature) quantum phase transitions? 

In this paper, we will investigate these questions by means of a low-energy approach where the 
symmetry breaking scheme is described by a Gk WZNW model perturbed by marginal current- 
current interactions with H invariance. The main result of this work is the fact that this general class 
of models possesses hidden non-trivial duality symmetries which are associated with the underlying 
competing orders. These dualities are emergent, and exact, in the continuum limit. They can be 
viewed as the naive generalization for continuous symmetry groups of the well-known Kramers- 
Wannier (KW) duality symmetry of the ID quantum Ising model. Some duality symmetries are 
not new and have already been found in specific models such as the two-leg electronic spin- 1/2 
ladders at half-filling [14, 15] or away from half-filling [16-18], in two-leg spin ladder with four- 
spin exchange interactions [19, 20], and also in ID multicomponent cold fermions [12, 13]. Most of 
them have been revealed within a one-loop renormalization group (RG) approach of the interactions 
with the emergence of a dynamical symmetry enlargement (DSE) [14]. Such DSE corresponds to 
a situation where a Hamiltonian is attracted under RG flow to a manifold possessing a higher 
symmetry than the original bare theory. The main interest of this DSE scenario stems from the 
fact that the isotropic RG ray is usually described by an integrable field theory. The integrability, 
in turn, leads to the description of the low-lying excitations and the determination of the physical 
properties of the corresponding spin-gapped phases [14, 21-23]. However, as it will be shown 



1 This definition of the group G as the microscopic degeneracy group of the non-interacting theory, G mlcro , is not 
fully correct. A more rigorous definition of the group G will be given in section II: it is the maximal symmetry 
group of the interacting sector of the theory, which is of course closely related to G mlcro . 
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in the following, the duality symmetries are not restricted to a two-leg ladder or specific to the 
one-loop RG calculation but correspond to a general powerful non-perturbative approach in ID. 
In particular, they can be classified and depend only on the algebraic properties of G and of the 
actual symmetry H of the lattice model. As it will be shown, the resulting list of dualities turns 
out to be in correspondence with the classification of Z2 graduations of semi-simple Lie algebra. 
The interplay between these dualities and the DSE phenomenon will also be clarified in this work. 
In simple words, the dualities correspond to the different possible DSEs of the problem. The 
existence of these dualities symmetries enables one to relate different possible competing orders and 
to shed light on the global structure of the phase diagram of weakly-interacting fermionic models. 
Of course, our work does not solve the general problem of classification of zero-temperature ID 
gapped phases but, at least, the duality approach provides a classification of a subset of gapped 
phases from symmetry arguments. Moreover, the nature of the quantum phase transitions between 
these orders can also be determined, within our approach, by investigating the physics along the 
self-dual manifolds. 

The rest of the paper is organized as follows. The general low-energy effective theory, that 
we will study in this work, is presented in section II. Section III contains our most important 
results with the definition and characterization of the duality symmetries. These hidden discrete 
symmetries are introduced from two different complementary approaches. A first definition stems 
from the covariance of the low-energy effective Hamiltonian. The characterization and classification 
of these dualities are then deduced in section III A. A second definition is presented in section III 
B in light of the DSE phenomenon where duality symmetries are viewed as exact isometries of 
the RG beta-function. The nature of quantum phase transitions between different orders is then 
investigated. In section IV, the duality approach is applied to weakly-interacting multicomponent 
fermionic models for incommensurate filling. The special example of a generalized two-leg ladder 
is then considered in section V. Finally, we conclude in section VI and some technical details on 
the duality approach are presented in an appendix. 

II. GENERAL LOW-ENERGY EFFECTIVE HAMILTONIAN 

Let us consider n chains of weakly-interacting lattice spin-1/2 fermions Cj iQ a = 1, .., N where 
N = 2n is the total number of internal degrees of freedom of the problem. The generic model, that 
we have in mind in this paper, is described by a lattice Hamiltonian of the Hubbard type with 
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contact interactions: 



H = "* E E (4aC l+ l, a + H.C.) + E E ^K) ^ai^.as'W W 



where t is the hopping term and the independent couplings (on-site couplings U{ a .y) depend on 
the precise model under consideration, that can be specified by fixing its physical symmetry group 
H. In this respect, many other terms could be added to model (1), without affecting the following 
discussion which is mostly based on symmetry considerations. For example, more general short- 
range interaction terms could be considered as long as they respect the symmetry H of the model. 
More concrete examples will be discussed in section V. 

In the continuum limit, the non- interacting spectrum around the two Fermi points ±k F is 
linearized and gives rise to left- and right-moving Dirac fermions ^ aL ( R ) [6, 7]. The left (rights- 
moving fermions are holomorphic (antiholomorphic) fields of the complex coordinate z = v f t + ix 
(t being the imaginary time and v F is the Fermi velocity): ^ aL (z) , ^ aR (z) . They are normalized 
by the following operator product expansion (OPE): 

*L(*)*te(*) - - 



2ir (z — uj) 

*L.(*)<MW) ~ - , 6ab y (2) 

2-7T [Z — LO) 

In the weak-coupling regime, the effect of the interactions is described by the following general 
low-energy effective Hamiltonian density: 

H = H + H int = ~iv F (*LA*aR - *LA* l) + 9a J£ d a AB J L B , (3) 

where velocity renormalization terms and interactions with the same chirality will be neglected 
in this work and g a accounts for the coupling constants of the fermionic models. In the follow- 
ing, repeated indices are summed throughout this paper unless stated otherwise. The continuous 
symmetry of the non-interacting Hamiltonian of Eq. (1) is U(iV) = U(l) x SU(iV) where U(l) is 
the usual global charge symmetry (cj ja — ► e^Cj ja ) and the SU(iV) group describes the continuous 
symmetry of the remaining degrees of freedom, dubbed in the following "spin" sector for simplicity. 
This basis, which singles out the charge degrees of freedom, is natural for incommensurate filling 
since a spin-charge separation is expected in the continuum limit so that model (3) decomposes into 
two commuting pieces: H = 7i c + 7i s . In this case, the charge degrees of freedom display metallic 
properties in the Luttinger liquid universality class [5-7]. The Hamiltonian in the spin sector de- 
pends only on the currents <J^ R j = T A (A = 1, .., iV 2 — 1) being the generators 



of SU(iV) in the fundamental representation (we choose the nor malization Tr (T A T B ) = 5 AB ) and 
reads as follows: 

n s = AqT { V N F + 1} ( j£ J A + J L A J L A ) + 9a J A d% B J B , (4) 

where, in the following, normal ordering is always implied for the non-interacting term. Model (4) 
is nothing but the SU(iV)i WZNW model perturbed by a current-current interaction. 

In contrast, there is no spin-charge separation in the half-filled case and the charge degrees of 
freedom cannot be disentangled from the spin ones due to an umklapp process. A good starting 
point is then to consider the maximal global continuous symmetry of the non-interacting lattice 
Hamiltonian (1) i.e. G = SO(2iV). Physically, the appearance of the group SO(2iV) can be 
motivated as follows: due to particle-hole symmetry present at half-filling, the 2N local states 
{c\ a \GS) ,c ia \GS) | can play a symmetric role, where \GS) is a particle-hole symmetric many- 
body ground state. In the continuum limit, the SO(2iV) symmetry can be revealed by intro- 
ducing 2N real (Majorana) fermions from the N Dirac ones: VP a L(R) = (£a + i £ci+./v)l(r) 

/V2. 

The non-interacting Hamiltonian of model (3) then takes a manifestly SO(2iV) invariant form: 
Ho = —iv F (^audx^an — iahdxiah) /2. Within this representation, the currents J^L^ in Eq. (3) 
express as fermionic bilinears: = *27r£ oL ( R )£& L (R)> A = (a, b) , 1 < a < b < 2N. 

The Hamiltonian (3) contains only marginal interactions - this is the main, though important, 
restriction to our study - and includes several competing orders which are encoded in the physical 
symmetry of the problem, the group H. In Eq. (3), these symmetries are taken into account by the 
set of symmetric matrices D = {d a } that commute with all elements of H: [cP,H] = 0. On top of 
continuous symmetries, the Hamiltonian (3) has discrete symmetries such as lattice symmetries, 
chiral invariance, charge conjugation and parity. In particular, those discrete symmetries are 
responsible for the matrices d a to be real and symmetric. In addition, model (3) should be stable 
and renormalizable under the RG approach and we thus require the set of matrices d a to be closed 
under anticommutation: {D,D} C D [24, 25]. 

In the following, we will generalize the problem and make abstraction of the underlying fermions 
to consider the G/% WZNW model perturbed by current-current interaction with Hamiltonian den- 
sity: 

H = H + H in t = 47r ( fc F + gV ) (^r^r + ^l A ^l A ) + 9a J A d a AB J B , (5) 

where A, B = 1, . . . , dim(G) and g v is the dual Coxeter number of G (for instance g v = N for G 
= SU(iV) and g v = 2N — 2 for G = SO(2JV)) [2]. In this paper, G is a regular simple Lie group 
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and for application to weakly interacting fermionic models (1), one has G = SO(2N) with the KM 
level k = 1 for model (3) at half-filling and G = SU(iV) and k = 1 for the spin degrees of freedom 
of model (4) away from half-filling. The WZNW description (5) has to be supplemented by the 
defining G^ KM current algebra: 

#W(«0 - 7^72 + ^ ABC — • ( 6 ) 
(z — w) z z — w 

where a similar relation for right-moving currents holds with the replacement z,w — > z,w. In Eq. 
(6), we choose a basis of the Lie algebra q of G for which the Killing form is the identity and f ABC 
denote the structure constants of g which are normalized according to: f ABC f ABD = 2g v 5 CD . 

In absence of interactions, i.e. when g a = 0, the theory (5) is confer mally invariant and 
its critical properties are described by the G^ WZNW CFT with central charge c(G, k) = 
k dim(G)/(/c + g y ). As far as global symmetries are concerned, the non-interacting model dis- 
plays a G L x G R symmetry, which is generated by the charges Q^e K \ = / dxjA^y One obvious 
effect of the interactions g a ^ 0, on top of reducing the symmetry down to H, is to couple the left 
and right sectors of the CFT: the continuous part of the physical symmetry group H of model (5) 
is generated by a subset of the diagonal charges (i.e. that rotate simultaneously the left and right 
sectors of the theory) Q A + Q A , so that more rigorously H = (H L x H R ) di . For generic g a ^ 0, 
we will investigate situations where the interaction is marginally relevant so that a spectral gap 
opens and the conformal symmetry is lost. 2 One important point about the interacting part of 
the Hamiltonian (5) is that it contains a special, G-symmetric ray, i.e. there exist some couplings 
g& such that g%d a = 1 (1 being the dim(G) x dim(G) identity matrix). If this were not the case, 
it would mean that the Hamiltonian could be broken into smaller, commuting pieces, built on 
subalgebras of G, to which our approach should be applied separately. Along the G-symmetric ray, 
the Hamiltonian (5) simplifies as follows: 

^ G = 4 , ( ;; gV) + tftf) + 9 j*Jt (7) 

which is the Gross-Neveu (GN) [28] or Thirring model [29] built on the group G at level k. On 
this special ray, the global symmetry of model (7) is given more rigorously by G = (G L x G R ) diag , 
which is generated by the diagonal charges Q A + Q A . This last remark finally provides us with 



2 An exception corresponds to the case where the marginal current-current of model (5) is described in terms of 
currents with different KM levels. Even though the perturbation is marginally relevant the conformal invariance 
is restored at an intermediate infrared fixed point. The so-called chirally stabilized spin-liquids are examples of 
this class of models [26, 27]. 
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the correct way to introduce rigorously the group G, namely the maximal symmetry group of the 
interacting model (5). 



III. EMERGENT DUALITIES 



In this section, we study duality transformations on the general model (5). After a brief survey 
of basic facts about the famous KW duality in the quantum Ising model, we start by introducing 
duality symmetries that map the set of theories defined by Eq. (5) onto itself, from a rather 
algebraic point of view. This definition will enable us to obtain a classification of dualities depending 
on the nature of the group G. Then, we study the interplay of such dualities with the phenomenon 
of DSE. For readers not interested in the technical details, a brief summary of our main results is 
presented at the end of this section. 



A. Dualities and covariance of the Hamiltonian 



Dualities are precious tools that help our understanding of strongly correlated systems. Given a 
set of theories depending on some parameters {g a }, with Hamiltonian H ga [(f)] and fluctuating fields 
{(p(x)}, a duality Q is a symmetry operation that relates different points in this set of theories, 
Hg a W\ = Hg a W\i where parameters as well as fields are acted upon by the duality: £l(g a ) = g a 
and = with O 2 = 1. This ensures that the whole physical content of the theory Hg a can be 
deduced from that of the theory H ga : correlators are invariant, ( . . . (/)(x) • • • )~ = (• • • <A(^) • • ■) ga , 
where (. . .) is the quantum average defined by Hamiltonian H ga . 

1. Warming up: Kramers- Wannier duality 

The ID Ising model in a transverse magnetic field provides us with maybe the most famous 
example of such a duality, the so-called KW duality. It is worth recalling basic facts about this 
model, since it reveals striking similarities with the more general situation studied later. Its lattice 
Hamiltonian is: 



H 



Yl + 97°i\ , (8) 



where af are Pauli matrices on site i. Defining the operators Fl((rf) = [x\ = Ylj^&j and f2(of 



fjbf = crfaf +1 allows a dual description of this model, H = ^ 



92^1^1+1+91^ 



. The KW duality 



ft simply exchanges the couplings, Q(g 12 ) = g 2 \- At the self-dual point g 1 = g 2 , the theory is 
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critical and described by a c = i CFT which is the free Majorana fermion theory [2, 6]. Moving 
away from the self-dual point amounts to adding a mass term to the fermions, yielding the following 
Hamiltonian density in the continuum limit: 

11) 

n = -- (£ a x £ B - + ^£r£l, (9) 

with m oc g 1 — g 2 - Then the KW duality exchanges the two massive phases of the Ising model: when 
m > 0, it is in its ordered phase, (o~ z ) ^ 0, (/i z ) = 0, while the disordered phase is characterized 
by 7^ and (o~ z ) = 0. On top of Majorana fermions, there are two other fluctuating fields 
in the theory, cr(x) and [J,(x), that are the continuum limits of the Ising spin a z , and of the 
disorder operator respectively. The KW duality has the following possible representation on 
the continuous fields: £ L — > — £ L , £ R — > £ R , a «-> /i which is indeed a symmetry of model (9) with 
m — > —m. Note also that if the original 7L 2 symmetry is broken (e.g. by adding a magnetic field 
along a z ), then the KW duality no longer holds. More precisely, if one then requires the model to 
be globally invariant under KW (up to a redefinition of the couplings), one is forced to include a 
new term in the Hamiltonian, namely a term proportional to /i 2 which is not local with respect to 
the original building blocks af. In other word, preserving the possibility of KW duality requires to 
enlarge drastically the class of perturbations. This is a general fact: as we will see, for the generic 
model (5), the set of allowed dualities strongly depends on the symmetry group of the model. 
Moreover, a quantum critical behavior is also likely to emerge at the self-dual points. In this 
respect, the KW duality approach has been fruitfully exploited in the past to obtain the massive 
phases and quantum phase transitions of two-leg spin ladders for weak interchain interactions [30] . 

2. Characterization and classification of dualities 

Motivated by the example of the Ising model, we now investigate whether the theories defined 
by (5) admit non-trivial transformations that relates different points in this set of theories. Since 
the models that are considered are built out of the current fields, it is natural to define such a 
transformation - call it Q, — by its action on the currents, J — ► J = Q(J'). In close parallel to 
the Ising model, we demand that the model defined in terms of the new variables J , corresponds 
to the original one up to a redefinition of the couplings, g — > g: 

H(g,J)=H(g,J). (10) 

We will see shortly that such transformations can exist, and, moreover, under very general assump- 
tions, are involutive; hence it is sensible to coin them dualities. 
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In this paper, we will restrict ourselves to the simplest case of transformations that act linearly 
on the currents. Imposing that the KM algebra (6) be invariant prevents left-right mixing of 
the currents, thus resulting in: <J^ R ) = ( u l(r)) AB 'J^ R y w l(r) being dim(G) x dim(G) matrices 
representing the action of on the currents. From the invariance of the OPE (6), w L ( R ) are necessary 
orthogonal matrices. By writing uj l = ujuj r , with u = uj^^ 1 , we recognize that the transformation 
can be decomposed into the product of a pure diagonal G-rotation lj r (i.e. that acts simultaneously 
on the right- and left- chiral sectors) and of a transformation acting only in the left sector, 10. It 
can be shown (see Appendix A 1 for details) that to and the diagonal rotation lo r have to leave the 
Hamiltonian globally invariant (i.e., they separately fulfill Eq. (10)). Therefore, the problem of 
finding all possible dualities factors out in two independent sub-problems. The diagonal rotation ui R 
corresponds to a mere change of basis in the representation of the Hamiltonian. If the low-energy 
theory that one considers is the continuum limit of some lattice model, this diagonal rotation is 
the continuum representation of a local, unitary transformation of the lattice operators. Here, we 
will not be interested in these diagonal transformations. We mention that such transformations 
have been studied on the lattice in the context of the two-leg Hubbard ladder [15] and in the 
two- leg spin ladder with ring-exchange interactions [19, 20]. One can convince oneself that they 
appear whenever there is some accidental degeneracy in the decomposition of the Lie algebra of G 
in irreducible representations of H. 3 

In the following, we will investigate the more interesting transformations that affect only one 
chirality sector, setting therefore lo r = 1. These transformations will be called dualities, and 
generically denoted by 0. The matrix uj is the representation of Q on the KM current fields; Q 
twists the current algebra: 

n-A H 7}A _ , AB rrB 

Jl * J L — U J L , 

JA _0, jA = jA_ 

Contrarily to the diagonal rotations, these transformations that affect differently the left- and 
right-sector of the theory cannot be viewed as the continuum limit of a local transformation of 
the underlying lattice fermions. In this sense, they are non-local, like the KW duality in the 
Ising model which involves non-local transformation of lattice operators. Moreover, in general and 
contrarily to the KW duality, there is no known, simple lattice realization of the dualities in terms 

li An archetypical example of this is the spin-orbital model, in which spin and orbital degrees of freedom are treated 
on an equal footing, both transforming under SU(2) rotations, with H=SU(2) or j :) xSU(2) S p m C G=SU(4). Then, 
there is a SU(4) automorphism exchanging spin and orbital degrees of freedom. We will come back to this example 
in section V. 
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of fractional variables (akin to the Ising spin and disorder operators in the case of KW duality, 
that are fractional with respect to the fermion) these transformations are thus in general emergent 
in the continuum limit 4 . 

We now proceed to characterize the set Ql of possible dualities O. It turns out that the structure 
of $) is remarkably simply extractible from few basis data: the physical symmetry group H, and 
the maximal symmetry group G of the interaction. 

First, one observes that uj should preserve the KM current algebra (6); it results that f2 is 
an automorphism of q. The covariance of the interaction term imposes further constraints: there 
should exist couplings g a = Q(g a ), such that 

g a d a = g a d a uj. (12) 

Multiplying this relation by constant symmetric matrices and taking a trace, we see that the 
transformed couplings g a are linearly related to initial couplings g a . Moreover, condition (12) is 
strong enough to completely characterize allowed dualities: it is possible to show (see Appendix 
A 1) that the set of allowed dualities is given by: 

^ = C(H)| inv , (13) 

i.e., the matrix to should be the representation (on the adjoint representation of G) of those elements 
of the center of H 5 , that are involutive, lo 2 = 1. This characterization of Q has direct, simple 
consequences. First, it shows that there is a finite number of dualities. Si is indeed isomorphic to 
the finite group (^2)™, n being some positive integer, and has 2 n elements. Furthermore, to each 
duality, there corresponds a symmetry operation S that belongs to the physical symmetry group H, 
which has the following properties: S 2 = 11 (S is a Z2 symmetry), and S commutes with all other 
symmetry operations of H. Explicitly, a representation of S can be given as: S = to x to, where the 
two operands of the tensor product act on the left- and right- sector of the theory, respectively. A 
duality therefore corresponds to the "square root" of some exact involutive symmetry of the lattice 
Hamiltonian that commutes with all other symmetries. 

Let us denote by Q a the different dualities, *2> = {Q a }- Using the properties [w„, wj] = and 
lu 2 = 1, we deduce that there exists a basis for the currents that simultaneously diagonalizes the 
action of dualities on the currents, the eigenvalues being ±1. In the following we will work in such 



4 Note that our study does by no means exclude the possibility of exact, complicated, necessarily non-local lattice 
realizations of dualities in some specific examples. 

5 We recall that the center of a group H is defined from the elements which commute with all elements of H. 
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a basis: 

Q a (J A )=e A J A , e A = ±l. (14) 

This yields a more physical interpretation of dualities. Introducing Noether charge and current 
densities associated to the G-invariance, J A = J A + J A and J A = J A — J A (the subscripts 
correspond to space-time indices, xq = t and x\ = x, and in a G-invariant theory dtJ^ +d x J A = 0), 
they transform as: 

n a (j£) = j£, n a (J A ) = J A , e A = +l, 

n a (j A ) = j A , n a (j A ) = j A e A = -i. (15) 

In a good basis, a duality just amounts to the exchange of the role played by Noether charges and 
currents associated to G-invariance. Thus, it is tempting to view dualities as generalizations of the 
U(l) duality that exchanges the role played by electric and magnetic charges (this is known as "T- 
duality" in the context of string theory [31]), to a non-Abelian theory, with the further constraint 
that they have to be compatible with the original H-invariance of the microscopic theory. This 
will have important consequences for the labeling of states by quantum numbers. Indeed, if the 
eigenvalues of the charge operator Q A = J dxJ^(x) happen to be (approximate) good quantum 
numbers in a phase and if e A = —1, it follows that the (approximate) good quantum numbers in 
the dual phase £l a (^) will be current quantum numbers: quasiparticles carry generalized currents. 

We now make use of a mathematical result on Lie algebras that allows for a complete classifica- 
tion of dualities. Given O a £ @, using Eq. (14), it is possible to decompose the Lie algebra of G 
into two orthogonal subspaces: = 0y ©0" , 0y (fl" respectively) being generated by those elements 
of q with e A = +1 (e A = — 1 respectively). Mathematically, this decomposition can be associated 
to any involutive automorphism [lo in our case), which is called a Z2-grading of q [32]. Each 7L<i 
grading Q is characterized by an invariant subspace 0y , which can be shown to have the structure 
of a semi-simple Lie algebra, and is defined by £l(X) = X, VX £ 0ii. One has the orthogonal 
decomposition (with respect to the Killing form) = fly © p, with Q(X) = —X, \/X £ p, and the 
Lie structure of breaks up as follows under the action of 0: 

[0||,0||] =0j|, [0||,P] =P, [P,P] =0||- (16) 

Now, there exists a complete classification of Z2-gradings for simple Lie algebras [33]. This means 
that once G is given, irrespective of the physical symmetry group H, one knows what are the 
different possible types of dualities: they are of the form to = Uu)U~ l , where U belongs to G, and 
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Type 


G 


All 


AI 


SU(iV) 


so{N) 


An 


SU(2iV) 


sp(2N) 


AM 


SV(p + q) 


S (u(p)®u(q)) 


BDI 


SO(p + q) 


so(p) ffiso(g) 


BM 


SO(27V) 


u(N) 


CI 


Sp(2iV) 


u(N) 


CI 


Sp(2p + 2q) 


sp(2p) sp(2q) 



TABLE I: Exhaustive list of Z2-gradings for simple Lie algebras. For the Z^-grading of the AW type, we 
use a somewhat non-standard notation to designate the algebra of block diagonal matrices, with blocks A 
and B on the diagonal, where A (respectively B) is a p x p (respectively q x q) Hcrmitian matrix with the 
constraint Tr(yl + B) = 0. In this case, is isomorphic to su(p) © s\i(q) © u(l). 

where cD is one representative of the possible, classified, Z2-gradings of q. For regular simple Lie 
groups G, the Z2-gradings are presented in Table I [33]. 

In section IV, we will investigate in detail the case G=SU(iV) for applications to fermionic 
models (3). In this picture, the role of the physical symmetry group H is to select what are the 
realized dualities amongst the set of possible ones. That the set 3l of dualities is not reduced to 
the trivial identity element will be clear later on the example presented in section V. 

3. Transformation of fields 

So far, one has fully characterized the set of dualities @, that depends only on the maximal 
symmetry group G and the physical symmetry group H. This was done by investigating the action 
of dualities on particular fields, the left-moving currents of the theory (the action on the right 
currents being trivial). To investigate the effect of dualities on observables, we need more generally 
their action on any local field of the G& WZNW model (the unperturbed theory); this action turns 
out to be particularly simple. 

The whole field content of the WZNW model can be recovered by considering the KM primary 
operators ^a l ,a b [2]. Since the WZNW model is invariant under independent global chiral G- 
rotations in the right and the left sectors, each primary field is labeled by two G-highest weights 
A L and A R that dictate how the primary field transforms under the group G L xG R . Each primary 
field is thus a tensor ($A L ,A R ) a a , «l(r) = 1) •••> dim(A L ( R )). The precise set of representations A L ( R ) 
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that occurs in the WZNW model defines the field content that depends on the level k [2]. 

Before we give the action of dualities on these primary fields, we need a few basic facts about 
the group of automorphisms Aut(G). It contains the subgroup of inner automorphisms Auto(G), 
that can be identified with conjugation by elements of G, i.e. those transformations A\j that 
act on G as: g — ► Ajj(g) = UgU -1 , with U £ G. Thus, inner automorphisms can be viewed 
as the group of changes of basis in the representations of G 6 . In fact, inner automorphisms do 
not exhaust Aut(G). Indeed, Aut(G)/Auto(G) is in general a non-trivial discrete group. When- 
ever it is not trivial, we will denote its non-trivial elements generically by r. This leads us to 
introduce outer automorphisms, that can be written as A o r with A EAuto(G), and to write 
Aut(G)={Aut (G),Aut (G)or}. Contrarily to inner automorphisms, that leave globally invariant 
any G-representation, outer automorphisms generically exchange representations, i.e. r(A) 7^ A. 
Outer automorphisms are in one-to-one correspondence (up to conjugation) with the symmetries 
of the Dynkin diagram associated to G [32]. 

For regular simple Lie algebras, involutive outer automorphisms do exist and will be relevant to 
our study (they will give rise to "outer" dualities). An exhaustive list of outer automorphisms is: 
for G=SU(iV), the non-trivial element r corresponds to charge conjugation, that maps a represen- 
tation on its charge conjugated partner - the corresponding outer Z2-gradings are of type AI and 
AF (see Table I) - ; for G=SO(2iV), the non-trivial element r corresponds to "SO(2N) parity", 
that exchanges the two spinorial representations - the corresponding outer Z2-gradings is of type 
BDI when p and q are odd (see Table I) -. 

Let us now describe the action of dualities on fields. Since a duality Q E @ affects only the left 
sector, its action on the primary field ^a l ,a r is readily found: it is simply given by the action of 
the automorphism uj on the left representation A L . If Q belongs to the inner class, it just amounts 
to a change of basis in each representation A L . Denoting by U\ L the matrix performing the change 
of basis, the primary operator becomes: 

$a l ,a r $a l ,a r = U Xh • $a l ,a b , inner. (17) 

In this expression, a matrix product is implied; in components it reads as follows: (^ ) A L ,A R ) a a = 
(f^A L ) OL b L (^A L ,A R )b L a • O n the other hand, outer dualities map in general the representation A L on 
another representation A£ = t(A l ), and also perform a change of basis in A*. Thus, the primary 



6 Strictly speaking, Auto(G) contains more than conjugation by elements of G, since it is the connected component 
of the identity in Aut(G). This subtlety will play no role in our study. 
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operator transforms under duality according to: 



$A L ,A, 



Q 



$A L ,A B = U X; ■ $A L *,A 



outer. 



(18) 



« ) 



Let us now briefly restate the main results of the preceding sections: 

• To each family of models (5) defined by (i) an unperturbed critical model, namely a WZNW 
model built on the symmetry group G, at level k (for the fermionic models (3), the level is 
fixed to k = 1) and (ii) its bare symmetry group H, there is a finite set of exact dualities Ql. 

• To each model belonging to this family, defined by its bare couplings g a , one can associate 
a finite number of models defined by couplings £l(g a ) where ft G 3> is an acceptable 
duality. These dualities are exact in the continuum limit. The physical properties of 
can in principle entirely be deduced from those of j% . For example, suppose the ground 
state of model ^# is specified by some order parameter O, that can be expressed in terms 
of a primary operator: O = Tr (MQ\ h where M is some matrix encoding the "texture" 
of the phase of model j& '. Then, the ground state for model is immediately given, with 



It follows that it would be sufficient to understand the physics of one "fundamental chamber" 
- a subspace of the space of the couplings g a - from which physical properties for the whole space 
of couplings could be deduced by acting with dualities. Nevertheless, understanding the physics of 
the "fundamental chamber" is still of course a complicate task. In particular, there is no reason 
for this "fundamental chamber" to display a single type of physical behavior - it may contain 
many different phases. In section IV, we will give a physical picture of part of the "fundamental 
chamber" in the SU(iV) case to exhaust its physical properties. 



In this section, we present an alternative characterization of the dualities found in the previous 
section: they are in one-to-one correspondence with the different ways the bare physical symmetry 
group H of the model can be promoted to a larger dynamical symmetry group G. This notion 
of DSE phenomenon, which can be made explicit by means of the RG approach [34], is recalled 
for completeness and applied to model (5). Then, important consequences for the global phase 
diagram are drawn, that allow for a general picture thereof. 




B. Dualities and DSE 
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1. DSE phenomenon 

Perturbing a ID critical model with a continuous symmetry offers in general a whole family of 
possibilities, that one can label by the symmetry breaking pattern and by the type of perturbing 
operator; in our case, once chosen the general class of perturbations of the current-current type 
(5), the model is specified by the pattern G~> H. Generically, a marginal relevant perturbation 
results in the opening of a spectral gap in every sector of the theory. However, in one dimension, a 
simplification can occur that might allow for a kind of universal description of these fully gapped 
phases. As we have seen in section II, the family of models (5) contains a special, isotropic ray 
(with couplings proportional to where model (7) displays an extended global G symmetry. 
This GN model, that we will call ^#o> has the following Hamiltonian density: 

n G = H + gJ*J L A . (19) 

For g > 0, this model is a massive integrable field theory and its spectrum is known for k = 1 
[35-37] and in the general k case [38-40]. The coupling g flows to strong coupling under RG, and 
a mass scale is dynamically generated: m oc e~ cst ^ 9 . 

It turns out that physical properties of the generic model (5) can be weakly sensitive to 
anisotropy, i.e. to departure from the isotropic ray, in a sense we are going to discuss now. To re- 
state it differently, the effective theory describing the bare anisotropic model at low energy, displays 
approximately a larger symmetry G, and the physics is that of the maximally symmetric model 
^#o- This phenomenon has been coined DSE. An indication for DSE can be gained from consid- 
ering the RG beta- function for generic models of the form (5). As shown in the detailed study of 
Refs. 34 and 41, the beta- function generically indicates that the more symmetric ray is attractive 
under RG (writing g a = gg^ + Sg a the quantity 5g a /g flows to zero at one loop lim.£^o(5g a /g) = 0, 
with I = ln(Aao/v F ) the RG time, A the running cut-off and ao the lattice spacing). The model 
thus flows to more and more symmetric theories at low energies. The RG flow defines a scale A DSE , 
akin to Josephson length, where anisotropy becomes negligible: (5g a / g)(ln(A DSE ao / '%)) -C 1. Now 
focusing on the low-energy spectrum (precisely speaking at energies E < E\ ow <C A DSE ), we can 
have a non-trivial regime described by field theory only if one considers the triple limit 

(strong DSE) m < E {ow < A DSE < v f (Iq \ (20) 

where m oc e~ cst ^ 9 is the typical gap scale of the system (in this asymptotic DSE regime a single 
coupling g is enough to describe the flow). Regime (20) is universal: in that limit the lattice model 
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(1) and model (19) share the same quantitative physical low energy (E < E\ ow ) properties (this 
quantitative DSE phenomenon could be checked explicitly in some cases where the anisotropic 
model itself is integrable, allowing for a complete determination of all physical quantities [34]). It 
is not very surprising that when dealing with gapped systems, the notion of universality requires 
far more restrictive conditions - as far as the regime of bare parameters of the underlying lattice 
model is concerned - than for critical systems (for which a Josephson length Aj can also be defined 
and that only require Ei ow <C Aj <C v F aQ 1 ). Note that (20) requires in particular m <C A DSE , which 
implies that the running coupling constants are small at the DSE scale: <7 a (A DSE ) <C 1. This shows 
that strong DSE is essentially a weak-coupling phenomenon [34, 41], and this justifies the one-loop 
approximation. 

It is also possible to define DSE in a weaker sense, in which one does not require that the 
physical quantities be quantitatively the same as those of the isotropic point. The meaning of 
DSE in the weak sense, can maybe be best grasped in light of adiabatic continuity: it simply 
reflects the fact that the ground state at the isotropic point can be adiabatically connected to the 
ground state of the (weakly) anisotropic models; moreover, this adiabatic continuity also holds for 
the quasiparticles that can still be labelled by G quantum numbers: each of the G-multiplets A 
of quasiparticles of the isotropic GN model splits under the (small) residual anisotropy breaking, 
m\ — > {m\, m\, ...), with a splitting 5m\ = max(|m^ — mjj) <C m\. This gives the hint that if one 
does not insist on the physical quantities to be asymptotically identical to those on the isotropic 
ray, then one can allow for finite, small deviations from the limit m/A DSE — > 0. This has been 
shown in the particular case of the SU(4) Hubbard chain at half-filling [42], where continuity could 
be checked numerically at finite anisotropy in a non-integrable model. Thus, DSE in the weak 
sense simply means that the group G is an approximate symmetry of the low-energy sector of the 
theory, that allows for a description of the low-energy sector - very much in the way approximate 
symmetries allow for a classification of baryons and mesons in particle physics. As shown in the 
following section, the splitting of the G-multiplet under the residual anisotropy can be controlled in 
a perturbative expansion, which means in particular that the quasiparticles of the isotropic point 
(of model ^#o) survive anisotropy 7 , and establishes the aforementioned adiabatic continuity on a 
firm ground. This weak-DSE regime, corresponding to: 

(weak DSE) m < A DSE (21) 



Of course, this cannot be the case irrespective of the perturbation type: our statement holds provided the anisotropy 
is due to current-current terms. 
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is characterized by adiabatic continuity of the low-energy spectrum (in particular, it is a smooth 
crossover and the symmetry of the ground state cannot change). Of course, symmetry enlargement 
cannot occur everywhere in the phase diagram: it has been established on the level of a perturbative 
expansion around the isotropic ray, and if the bare anisotropy is too strong, one expects it to break 
down. In particular, when A DSE /m decreases and Sg/g becomes of order one, level crossings start 
to occur in the lowest energy sector of the theory and the labeling in terms of G quantum numbers 
losses its meaning. In the following section, we will actually see that the existence of dualities, that 
where introduced in section III A purely on group theoretical considerations, implies that is has to 
break down, and in a rather elegant way, since it will allow to grasp at a glance the (very rough) 
structure of the phase diagram. 



2. Interplay with dualities 

We are now in position to make the connection with the dualities introduced in the previous 
section. It is possible to sketch the phenomenon of DSE as follows: 

Gl x Gr interaction , H = (H L x H R ) diag DSE , G ffi , (22) 

where the global symmetry group is indicated (from left to right, respectively) before perturbation, 
in the bare H-symmetric model (5), and in the effective low-energy theory in the infrared (IR) 
regime after DSE has occurred for instance in the weaker sense with the symmetry group G IR . 
One can then ask the following question: is diagram (22) unique? 

It turns out that the answer is positive up to twists in the current algebra of the form (11), i.e., up 
to dualities. Indeed, the final point of diagram (22) is in general of the form G IR = (G L x Gr)^ 
with generators given by J dx {j A + J^) , where J A = uj ab ' J b and u G 3> is a duality (see 
Appendix A 2 for a proof). Dualities can thus be viewed as the different possible DSEs compatible 
with a given bare symmetry group H. To each duality Q E @, there corresponds a generalized 
isotropic, G-symmetric ray, described by the following Hamiltonian density: 

Hn = H + gJ*LU AB J L B , (23) 

a model that we will denote Being related by a duality to the fundamental isotropic model 

(19), it is also an integrable field theory. 
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3. Symmetry enlarged phases 

Since dualities are exact isometries of the beta-function, model (23) is also attractive in the RG 
sense, so that DSE also occurs in the vicinity of the (generalized) isotropic ray. Now, the fact that 
DSE breaks down if the bare anisotropy is too large acquires a new meaning: it has to be so, since 
the different DSE schemes corresponding to different dualities are incompatible; in other words, 
the different symmetry enlarged phases compete. We will turn shortly to the study of this 
competition. 

What about the robustness of DSE, i.e. what happens when one moves away, in the space 
of bare couplings, from the isotropic ray = ^(Sa) (with g^d a = u, see Eq. (12))? One way 
to answer this question is to consider the one- loop RG flow, with initial conditions g a (a-o 1 ) = 
9( a o 1 )(da + $a{ a o $a <C 1. The latter has to be cut-off before the couplings become of order 
one, i.e. at a scale A > m, where m is the mass scale. Departure from the isotropic ray then 
results in a residual anisotropy, the couplings being g a (A) = g(A)(g^ + S a (A)). The study of Ref. 
34 ensures that 5 a (A) scales to in the scaling limit Aao << 1 when approaching the isotropic 
ray, resulting in corrections of order (g(a.Q 1 )5 a (aQ *)) in the physical quantities, with j3 > an 
exponent that depends on the groups G and H under consideration. 

An alternative argument, that establishes the adiabatic continuity of the low-energy physics 
when departing from the isotropic rays, can be given using on the integrability of the isotropic 
model (23). At scale A, the model can be described by the following Hamiltonian density: 

H = Tin + g(A)5 a (A) J R A d a AB J L B , (24) 

where the last term is the residual anisotropy. It is possible to investigate, by form factor methods, 
the effect of this term which constitutes a perturbation on top of an integrable model [10]: if 
the perturbation is local with respect to the quasiparticles, then the quasiparticles survive the 
perturbation, i.e. there is adiabatic continuity in the spectrum [10]. Since the perturbation is built 
on currents, that are the local generators of conserved charges, one concludes to mutual locality 
between the perturbation and the quasiparticles. It results that the low-energy spectrum of model 
(24) is described by the quasiparticles of the integrable model (23). 

4- Quantum phase transitions 

Let us now consider two different symmetry enlarged phases. The low-energy physics of each 
of them is captured by generalized isotropic models of the form (23), to be called and ^#q 2 . 
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These two models are exchanged by the duality Q = (recall that dualities form a group). 

In the space of bare couplings, amongst the different trajectories connecting these two isotropic 
rays, there is a special, maximally symmetric one, that allows to study, in a minimal model, the 
competition between the two phases and ^n 2 . Again, we invoke DSE to argue that this 

maximally symmetric trajectory is attractive under RG. 

Recalling that the duality $7 corresponds to a Z2 grading of the Lie algebra of G (g = gn © gj_), 
there exists a basis \jJn% J\ } for the currents with the following properties: 

^(^|l(r)) = ^\\l(r) ^(^J_l(r)) = T^IlCr), ( 25 ) 

where a = 1, . . . , dim(g|| ) and 6=1,..., dim(gj_). Then, the maximally symmetric model interpo- 
lating between the phases an d ^h 2 has symmetry G11 and is given by 8 : 

Wx_ 2 = n + g\\ J^fj\\L + 9± JL^i'JL (26) 

When choosing = g±, one recovers model ^h x , whereas g\\ = —g± yields model ^h 2 . There is 
a special point where the competition between the orders is at its climax: the self-dual point ^*„ 2 
corresponding to g± = 0, or, in term of the original couplings, to Q(g a ) = 9a- 

The model (26) can be further simplified by acting with the duality fii: this amounts to study 
the competition between the phase ^\ = ^#0 an d ^h, and in (26) to replacement lo\ — > 1, which 
will be assumed in the following. Then, at the self-dual point, the Hamiltonian density takes a 
simpler form: 

nh = n + 9 j( R Ji, (27) 

i.e., the currents J7j_ do not appear in the interaction term. This means that there is room for 
criticality at the self-dual point, with the decoupling of some degrees of freedom. The currents 
generate an affine algebra gy at level k* = tk, where l, a positive integer, is the Dynkin index for the 
embedding gii C g [2]. It is then possible to separate the "parallel" and "perpendicular" degrees of 
freedom at the level of the non-interacting Hamiltonian: TCq = TCo(q, k) = 7Yo(fl|| > k*) + Ti-Ai where 
the last piece is the Hamiltonian of model A describing the degrees of freedom that decouple at 
the self-dual point. Model A (that can be trivial, i.e., with vanishing central charge), is the coset 
model 

A = -^, (28) 

8 In the case G|| is not a simple Lie group, there are in general several couplings g\\, a - This does not affect our 
analysis. 
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with central charge = c(G,k) — c(G\\,k*), where we recall that c(G,k) = k dim(G)/(fc + g v ) is 
the central charge of the WZNW model. Model A can be non-trivial, with non-integer central 
charge. Concrete examples in the context of fermionic models will be given in sections IV and 
V. Of course, is only an upper bound to the central charge of the "perpendicular" criticality 
at the self-dual point: one has to check that no other relevant operator in model A is allowed by 
symmetries, that would spoil criticality. One also has to check whether the parallel sector itself 
can develop criticality. 




FIG. 1: Sketch of (a 2D-projection of) the phase diagram for the example of a model with four dualities: 
Qq = 1, Oi, O2 , and = ilif^- Bold lines correspond to symmetry enlarged rays ^#n, dashed ones to self- 
dual models ^n-fi' (f° r clarity some of the self-dual models have been omitted). All lines meet at the G& 
WZNW critical point (i.e. non-interacting point). Gray areas indicate regions of DSE and the dotted lines 
bounding it marks a cross-over to a regime where DSE looses its meaning (i.e. the splitting of G-multiplets 
due to anisotropy becomes of the order of the gap scale to). Four representative points together with the 
dualities exchanging them are also depicted. The dashed curve 7 symbolizes the path in parameter space 
which corresponds to the minimal theory (26) interpolating between phases ^#1 and . See figure 2 for 
a qualitative picture of the low-energy spectrum along this path. 



We thus arrive to the following general picture of the phase diagram for model (5) - see figure 

1: 

• There is a finite set of generalized isotropic rays where the bare model displays an exact 
enlarged G-invariance. Each ray is labelled by a duality Q 6 @, and the global invar iance 
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group G is obtained from G = (G L x G R ) dia g by acting with 17 on the currents (see Eq. (11)). 
The Hamiltonian on those rays is given by Eq. (23). 

• Each of those isotropic rays is attractive in the RG sense, and to it there corresponds a 
pocket with finite extension in the space of bare couplings. In this pocket, the low-energy 
sector of the theory is adiabatically connected to that of the isotropic model. The low-energy 
spectrum can be described by making use of the approximate G-invariance. In particular, 
the integrability of the resulting model gives a description of the massive quasiparticles which 
are organized in G-multiplets. 

• When the bare anisotropy is too large, DSE cannot hold any longer, i.e. the energy splitting 
of the G-multiplets due to anisotropy becomes of the order of the gap scale m: the low-energy 
physics cannot any longer be described in terms of a single energy scale m, and the situation 
is in general complex. 

• In between the symmetry enlarged pockets lie self-dual manifolds, which are the loci where 
anisotropy is maximal. A simplification occurs, with in general a decoupling of degrees of 
freedom, accompanied with criticality. The quantum phase transition that results - which 
in the simplest cases describes the transition between the DSE phases - is captured by the 
minimal model with Hamiltonian (26). Figure 2 presents a qualitative picture of the low- 
energy spectrum of this minimal model, along a path joining two distinct symmetry enlarged 
phases. For clarity, one chooses a path in parameter space along which the mass spectrum 
scales in a convenient way: apart from possible phase transition points, the typical mass 
scale is held constant. 

IV. FERMIONIC MODELS 

In this section, we come back to the fermionic models (3) defined in terms of N chiral Dirac 
fermions ^ aL ( a ). We have shown, in the previous section, that, irrespective of the bare symmetry 
group H of the model, there exists a finite number of possible dualities that can be classified. To 
each of them there corresponds a symmetry enlarged phase, which we now characterize. There 
will result an exhaustive list of DSE phases possibly supported by models (3). For each of those 
phases, the duality symmetries together with the integrability yield precious informations on the 
low-energy physics. As explained above, the maximal symmetry group supported by those models 
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self-dual point 



FIG. 2: Sketch of the low-energy spectrum of the minimal model (26) interpolating between two symmetry 
enlarged phases (see the path 7 in figure 1). 



is G = SU(iV) or G = SO(2iV), depending on the filling. Here we analyze the general G = SU(iV) 
case which is relevant to weakly-interacting fermions away from half-filling. The special N = 4 
case, i.e. two-leg electronic ladders, will be described in detail in section V. The general half-filled 
case with G = SO(2iV) is much more complicate to analyze and will be investigated elsewhere. 

A. Different classes away from half-filling 

For incommensurate filling, there is no umklapp process coupling spin and charge degrees of 
freedom: TL = Tt c + TC s . The charge fluctuations are decoupled and described by a massless bosonic 
field $ c and its dual C field, with Ffamiltonian density: 

which belongs to the Luttinger liquid universality class [5-7]. In Eq. (29), v c and K c are the 
charge velocity and Luttinger parameter, respectively, that depend on the microscopic details of 
the Hamiltonian and should be regarded as phenomenological parameters. For a commensurate 
filling, like one-electron per site, additional umklapp processes may appear depending only on the 
charge degrees of freedom and the low-energy charge Hamiltonian (29) becomes a sine-Gordon 
model. A charge gap might open for a sufficiently strong value of the interaction which results in a 
Mott transition. The non-trivial physics of model (3) corresponds to the remaining "spin" degrees 
of freedom which are governed by model (4) with G = SU(iV). 

According to the analysis of section III, we know that in the spin sector, there can only be a 
finite number of symmetry enlarged phases, that are classified by the SU(iV) Z2 gradings listed in 
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Table I. Due to DSE, each of these phases has a finite extension in the space of bare couplings, 
and its low-energy sector is described by the representative Hamiltonian (23), which is labelled 
by an involutive automorphism O of SU(iV). This results in four classes of symmetry enlarged 
phases corresponding to the trivial automorphism f2 = 1, and to the three automorphisms 
belonging to the classes AI , AM , AM of Table I. Recall that those are only the possible dualities: 
for them to be realized in a given model defined by the symmetry breaking pattern SU(iV) H, 
H must contain in its center an involutive element that gives rise to the duality. As we will see, 
physical properties show a strong distinction between symmetry enlarged phases ^(q according to 
whether they correspond to inner or outer automorphisms (here AI and AM correspond to outer 
automorphisms, and AM is inner). While for inner dualities, the order is of the density- wave type, 
phases associated to outer dualities display off-diagonal order, that comes along with superfluidity 
in SU(iV) spin space. In addition, we will see that outer dualities conjugate, i.e. maps fermions 
onto holes. A representative for these three classes for the SU(A^) group, in a suitable basis for the 
fermions ^ aL (^ a R being invariant in the definition of the dualities (11)), is given by: 

• AI class : ^/ ah — > Jl(\l/ OL ) = ^j^. The invariant subspace is bo(N), with a basis defined by 
the subset of SU(iV) generators W aL T^ b ^> bh , T % being real and antisymmetric. Obviously, this auto- 
morphism conjugates representations, i.e. it maps the representation A onto its charge conjugated 
partner A*. The lattice symmetry (the element of C(H) Jnv ) giving rise to O is charge conjugation: 

• AM class : ^> aL — ► J7(^ / aL ) = J a b^/ bL . This duality exists only for even N and the matrix J (J 
being the Sp(iV) metric) is given by —ia 2 x ljv/2- The invariant subspace gy is sp(N), with a basis 
defined by the subset of SU(iV) generators W aL T^ b ^f b such that JT l J = t T % . This automorphism 
also conjugates representations. The N = 2 is very special since the duality corresponds to the 
identity for the currents: Sp(2) ~ SU(2). The lattice symmetry giving rise to O is an "antisymmetric 
charge conjugation": c ja -> Y.b J abc) b . 

• AM class : \I> aL — ► Q(^ aL ) = (Ip^)^ 1 ^^, where I M (0 < p < q < N) is the diagonal matrix 
with p entries 1 and q entries -1. The automorphism belongs to Auto(G), i.e. it is not an outer 
automorphism. One has Qu = s(u(p) © u(g)) as an invariant subspace. The lattice symmetry that 
gives rise to this duality is: c- a — > —c- a for a > p. Dualities in the AM class are thus associated 
to a symmetry breaking pattern G ~> H that splits the ./V fermionic chains into two bunches of p 
and q = N — p chains. 
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B. Physical properties of the fundamental phase 

The physical properties of the symmetry enlarged phases can be deduced from those of the GN 
model (19), the representative Hamiltonian corresponding to the phase M\ labelled by the trivial 
duality = 1. 

• Ground state and spectrum: Its exact spectrum is known from integrability [35] and consists 
in N — 1 branches of quasiparticles transforming in the Young tableau A r with r boxes and one 
column (r = 1...N — 1), with masses m r = msin(irr /N). The quasiparticles with label r ^ 1, N — 1 
can be seen as boundstates of the "fundamental" quasiparticles with r = 1 and r = N — 1. The 
states are labelled by SU(iV) quantum numbers, i.e. by the eigenvalues of the N — 1 diagonal 
Cartan generators + Q°, a = 1, . . . , N — 1. 

The order parameter of this spin-gapped phase expresses simply in terms of the SU(iV)i WZNW 
matrix field g which is defined by [4]: 

g a6 = e-V¥*c*t A; (30) 

which is a pure spin field since it commutes with all operators of the charge sector. This matrix 
field identifies with the "fundamental" primary operator, g = &\ N _ 1 ,\ 1 - The order parameter reads 
then as follows: 

(Oj) = <Tr(g)> ± 0. (31) 

An order parameter on the lattice is readily found by considering the 2k F oscillating part of the total 
density operator rii = X^Li c la c i a> wnose continuum limit is related to g: n,2k F ~ e"^~^^ c Tr(g). 
This operator develops quasi-long-range order with a power-law behavior for the equal-times two- 
point functions: 

(n\ k¥ {x)n 2k M) (32) 

If the charge sector develops a gap by some process (for instance a pure charge umklapp term 
at commensurate filling 1/N (one electron per site) [43]), the resulting phase is a Mott insulator 
of the charge-density wave (CDW) type with {ri2k F + H.c.) ^ 0. This order may coexist with a 
2k F bond-ordering or spin-Peierls phase which is defined as: Osp = ^2j :a e~ 2 * ¥ ^ a ° c j, a c j+i,a 

+ H.c. 

A lattice order parameter involving only the "spin" degrees of freedom can also be found by 
considering a generalized dimerization operator in SU(iV) space: On = Ylj a 5 , /<S , j+ie _2ifcFja °i 
where = c'j a T^c- b is the SU(A r ) spin lattice operator on site j. The latter operator has a 
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continuum limit involving different harmonics [43, 44]: 

N-l 

Sf ~ y e 2irk * x a r AT r A {x), (33) 

where a r are vertex operators depending only on the charge boson - if the charge sector is gapped, 
they can be replaced by non-universal amplitudes (a r ) - and M A = Tt(T A $\ N _ rt \ r ) for r ^ 0, 
T A being the generators of SU(iV) in the representation A r (dim(A r ) x dim(A r ) matrices). The 
uniform spin density (i.e. r = 0) is M A = J A + J A , with a universal amplitude ao = l/2ir. Each 
of the harmonics Af A transforms in the adjoint representation of SU(iV) under the global diagonal 
rotations. When computed, the continuum limit of Oo shows that the SU(iV) dimerization operator 
develops a non-zero average value: 

(O d ) oc (Tr(g)) ^ 0. (34) 

We thus conclude, from the spin sector point view, that phase corresponds to a generalized 
dimerization phase, with an TV-fold ground-state degeneracy for a commensurate filling of one 
electron per site, which is formed from local SU(iV) spin-singlet condensation. 

• Susceptibilities and rigidities: Due to the spin gap, the susceptibility at zero temperature 
vanishes in all directions in spin space. Indeed, coupling the model to some field h A fdx(J A +J A ), 
one has 

d 2 E[h] _ n 

Xab ~ ~dhAdh? ~ U ' [6b) 
where E[h] is the energy of the system with the field h. In contrast, the quantity dual to suscep- 
tibility, namely the rigidity that measures the response of the system to twists in the boundary 
conditions (Cj + £/ ao a = Mab c- b with W [7] = e l ^ ATA ) for a system on a ring of size L, is non- vanishing 
in all directions of spin space. The rigidity is defined as p AB = q aq^b an d one has paa 7^ for 
all A = 1, . . . , iV 2 — 1. This is easily seen by bosonizing model (19): first one chooses a Cartan 
basis (H a ), a = 1, . . . , N — 1 with the first element H 1 along the direction A, and one introduces 
N — 1 bosonic fields </? aL ( R ), a = 1, . . . , N — 1, whose gradients represent the currents in the Cartan 
directions: i7^? R ) = V4^rd x tp aL ( R y The non-interacting part of the Hamiltonian (19) is that of N — 1 
free bosons TLq s = u^2 a ((d x (j) a ) 2 + (d x 6 a ) 2 ) /2, with (j) a = ip aL + tp aR the total boson field and 
6a = Pah — faR its dual field, while the interacting part of model (19), Tl[ n i s , is built from a sum 
of products vertex operators Vg a s a = e %t3alfah+% ^ aiPaR . The twist U = e %lHl , along the direction 
A in the boundary conditions, results in the following boundary conditions for the bosonic fields: 
4> a {x + L) = 4> a (x) and 8 a (x + L) = a {x) + ("f /^/^T)S a^ l. Now, due to global (diagonal) SU(iV) 
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invariance, and in particular to the invariance under the transformations generated by the Cartan 
H a that are represented by a shift on the bosonic fields, (faUn) ~^ ^«l(r) ± C a , Wi n t jS can involve 
only vertex operators with (5 a = (3 a , i.e. it does not depend on the dual field 8 a . It results that 
the twist in the boundary condition can be absorbed via the following canonical transformation 
affecting the boson 1 only: 6\{x) = 6\{x) — x, 4>\{x) = 4>i(x), yielding 

Paa = ^ + 0, (36) 

which signals the emergence of a finite rigidity. 

• Spin excitations: Another quantity of interest characterizing the phase is the dynamical spin 



structure factor S(q,u) = f dxdt e~ %qx ~ %wt \Sf+ x / ao (t)Sf(0)J . Using the continuum description of 
the spin operators (33), the Lehmann representation of the zero-temperature spin structure factor 
involves the form factors 

MO = (o\K A \t) , (37) 

where |£) is an eigenstate of the GN model (19). At low frequency, the behavior of (37) is fixed 
by one-particle states |£ r /) transforming in the representation A r /. Since the SU(A^) singlet repre- 
sentation does not appear in the tensor product of the adjoint by A r ', we deduce a strong selection 
rule: J r r (Cr') = Vr, r'. The consequence is that the spin structure factor displays no sharp peak 
structure at the quasiparticle poles q ~ 2rk F and uj ~ m r . This situation is sometimes referred to 
as incoherence of the spin excitations. 



C. Effect of dualities 



Acting with a duality 0, on the GN model (19) yields the generalized GN model (23), which, 
in spite of its simple connection to the original GN model, displays qualitatively different physical 
features. We now deduce the non-trivial consequences of the duality for the SU(A r ) group on the 
low-energy physics of the symmetry enlarged phases. 

• Ground state: For inner dualities, using Eq. (17), we observe that the WZNW field g = 
3>Ajv_i,Ai is mapped onto itself g — ► U\ N1 g. As a result, the order parameter of this phase is of 
CDW type, with a texture depending on the precise form of the duality symmetry encoded in the 
change of basis matrix U\ N _ 1 : 

(On) = (Tr(U XN _ lg ))^0. (38) 
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In a metallic phase, the 2k F component of the corresponding lattice operator = 
Y,ab c ], a ( U ^N~i)baC jb develops quasi-long-range correlations: 

(ngMn%l(0))~x-™*/ N . (39) 

If some umklapp process opens a gap in the charge sector, the ground state displays a CDW 
ordering characterized by (^2fc F ) ^ ®- ^ ne lattice generalized SU(iV) dimerization operator is now 
Op = A e A S^S^ +1 e~ 2zkF i a o with e A = ±1 according to whether the direction A belongs to gy 
or not. Using the result (34), we now conclude that ~ Tr(U\ N _ g), i.e. it develops a non-zero 
expectation value indicating a bond ordering. 

In contrast, for outer dualities, using Eq. (18), we see that the WZNW field is now changed 
into its dual, g — > U^g where: 

g ab = eV n o c m 

As a result, the order parameter of this phase becomes: 

(On) = (Tr([/ Al g)) / 0. (41) 

This translates into non- diagonal ordering for the lattice fermions and the order parameter becomes 
of a pairing type. The uniform component of the pairing density p^ = Ylab c j a (^Ai)&a c j b develops 
critical correlations 

(p n (x)p n (0)) ~ x - 2 l NK °. (42) 

In fact, phases associated to outer dualities spontaneously break a discrete symmetry of the original 
lattice Hamiltonian (1): the center Z^v of SU(iV), that consists of a discrete gauge redefinition of 
the phase of the underlying fermions Cj a — > e 2m / N Cj^ a . In the continuum limit, this transformation 
corresponds to special phase factors on the Dirac fermions: 

%N- *„ L (H)^e 2i7r/JV *aL(R)- (43) 

For iV > 2, p^ and Oq are clearly not invariant under this discrete symmetry. The spontaneous 
breaking of the center group symmetry is thus responsible for the formation of the spin gap of this 
phase. Note that it is likewise sensible to interpret the formation of a gap in phases associated to 
inner dualities as the spontaneous breaking of a dual group Z^r: 

*aL( R )^e ±2i7r/Ar *a L (K), (44) 
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and this Zjv symmetry has no local representation on the original lattice fermions. 

• Spectrum: The energy spectrum of model (23) is still that of the GN model (19), but now 
the charges labeling states, associated with the global symmetry G IR = (G L x G R ) diag , are the 
N — 1 Cartan generators Hfi = Q° + Q", a = 1, . . . , TV — 1, with = e a Q®, e a = ±1 according 
to whether the direction a belongs to gy or not. It is possible to show that for inner dualities, one 
can always find a Cartan basis untouched by $7, i.e. with e a = 1 Va, and we will work in such a 
basis. It results that states can be labelled by Noether charges Jq, see Eq. (15). For outer dualities 
though, it is not the case (for otherwise any representation A would be mapped to itself). Instead, 
there is at least one direction a in Cartan space such that Hq = J", i.e., the Noether current. We 
conclude that in the diagram (22), there is some sort of transmutation from Noether charges to 
Noether currents when going from UV to IR. 

• Susceptibilities and rigidities: Following again the same line of reasoning that lead to the 
results (35,36), one sees that dualities, in spite of the presence of the mass gap in every sector of 
the spectrum, results in the existence of directions in su(N) with non-vanishing susceptibilities 
Xaai an d vanishing rigidity P^a, namely those "transmuted" directions for which e A = —1. This 
is obviously consistent with the fact the associated quantum number along the direction A is a 
current and not a charge. If the duality is outer, any Cartan basis (that is, any choice of the 
quantization axis for the quantum numbers labeling the states) will be touched by the duality, i.e. 
it will contain at least one direction with e A = —1, and vanishing rigidity in that direction. It 
turns out that this is connected to an observable physical effect, that allows to distinguish between 
inner and outer dual phases: as shown in Ref.12, the breaking of the Z^r symmetry (43) in outer 
phases leads to a "confinement" of the current, elementary gapless excitations 9 carrying currents 
quantized in units of Ne. For a finite system put on a rotating ring, this results in a modified 
periodicity of the groundstate energy as a function of tangential velocity [12, 45]. 

• Spin excitations: For inner dualities, the conclusions drawn in the case of the trivial duality 
= 1 are not affected: spin excitations are still incoherent. However, for outer dualities, the spin 
density M A at wave vector q = 2rk F in Eq. (33) is mapped onto S~l(J\f A ) = Tr(T^ U\ r $\ rj \ r ), and 
owing to the tensor product A r (g> A r = A2 r ® • • • (indices are understood modulo N), the form 
factor of Eq. (37) J- r (^N-2r) 1S non- vanishing. Hence, the dynamical spin structure factor displays 
a sharp peak at q = 2rk F and u = m r : spin excitations in a charge-gapped system are now coherent 



9 This effect appears at incommensurate filling, i.e. with no gap in the charge sector. 
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. Outer dualities thus have a simple spectroscopic signature. 



Duality fl & phase 


order parameter 


inn/out 


l^ttipp ^vmmptrv 

giving rise to f2 


Trivi 3 1 

111V 1C\L 

SU(N) CDW-SP 


r . P ~2ik F ja 
C 3,a C j,a e 

r \ r p— 2ifc r ja 
j,a i+l,a c 


inner 


identity 


AI 

Symmetric pairing 


C j,a C j+l,a 


outer 


c — ► 

3, a J, a 


Aff 

Antisymmetric pairing 




outer 


(j = -jo-2 ® ljV/2) 


AM 

"Solit chains" p ^ ; 

g = N-p { : 


\p L~da<p q Z-*/a>p ) j,a^j,a^ 

/l V 1 V W f p -2ik F ja 
\p L-,a<p q l^a>p ) < 'j,a t 'i+l,a c 


inner 


C j,a ~* ~ C 3,a 
(o>p) 



TABLE II: The different classes of possible symmetry enlarged phases supported by the TV-component 
degenerate fermions in the absence of spin-charge coupling (incommensurate case). For each duality f2, we 
give the corresponding lattice order parameter, the nature (inner or outer) of the duality, and the lattice 
symmetry giving rise to f2. 



In Table II, we present a summary of the different possible symmetry enlarged phases for the 
A-component degenerate fermions close to the SU(AQi critical point, considered away from half- 
filling. 

D. Cross-overs and quantum phase transitions 

In this section we investigate self-dual points and their vicinity. We have seen in section IIIB4 
that those are points where the competition between two different symmetry enlarged phases is 
maximal, and can possibly lead to a quantum phase transition. For each pair of dualities fli 
and r^2, the region lying between phases ,^#ni and ^#r> 2 associated to them is best described in 
terms of the simplest model (26). By acting with the duality , one can reduce the family of those 
interpolating models that need to be investigated: it is enough to consider the models interpolating 
between and ^#C2, where ^ 1 is any non-trivial duality: 

Wi-n = Ho + g\\ + g± j£ ± J* ± . (45) 

10 This conclusion only holds in a charge gapped phase. If charge excitations are gapless, then the fluctuations in 
the charge sector (because of the presence of the operators a r in the spin operator, see Eq.(33)) will wash out the 
coherence in the spin-spin correlation function [22, 46]. 
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At the self-dual point (g± = 0), the interaction in the parallel sector gives rise to a gap for gn > 
(this is the case we will consider in the following 11 ): the spectrum in this sector is this of the 
Grii t GN model, where i is the embedding index of the embedding Gn C G associated to ft. There 
are, however, other degrees of freedom in general, that are not affected by the interaction: this 
decoupled sector is described by the coset model A = ^q^ 1 (see Eq. (28)). If model A is non- 
trivial, i.e. with central charge ^ 0, this is an indication that there is room for a continuous 
phase transition, in the universality class of the CFT A. If it is trivial, a phase transition might 
occur as we will see. The inner or outer nature of the duality will have different striking effects. 

• AI transition: The Dynkin index for the embedding SO(iV) C SU(iV) is i = 2, so that model 
A is trivial: = 0. The only degrees of freedom at the self-dual point are thus those of the 
parallel sector, which are governed by the following Hamiltonian: 

H* n = H + \Jf h Jf a , (46) 

where arc SO(iV)2 currents. 12 This is the k = 2 SO(iV) GN model which is fully gapped for 
A > 0. 13 In particular, it has been shown in Ref. 47 that this model is related to the Andrei-Destri 
[49] integrable model which is exactly solvable by means of the Bethe ansatz approach [49, 50]. 
Therefore, no criticality emerges at the self-dual point. However, recalling that 0, is outer and 
conjugates representations, we can show that the order parameters on both sides cannot coexist: 
while Tr(g) is invariant under the center of SU(N), that is under the Z^v transformations (43), it 
is not the case for the off-diagonal order parameter Tr(L r A 1 g) of Eq. (41) for N > 2. The Zjv 
symmetry is thus spontaneously broken in phase ^q. The opposite is true for the dual symmetry 
Zjv (44), which leaves invariant g but not g, so that it is spontaneously broken in phase We 
conclude to the occurrence of a first-order transition. 



The case g\\ < 0, for which the interaction at the self-dual point is marginally irrelevant, will not be analyzed 
here. It corresponds to situations where additional criticality is present at the self-dual point: criticality of the 
multicritical point extends in parameter space. We will present in detail an example of this in section V. 

12 We consider here the generic case N > 2 since the N = 2 case is very special. As discussed in Ref. 47, the resulting 
transition for N = 2 is described by a SU(2) self-dual sine-Gordon model which displays an U(l) quantum criticality 
[48]. 

13 A subtlety arises here from the fact that the conformal embedding SO(iV)2 C SU(iV)i, yields a non-diagonal 
invariant for the partition function in the SO(iV) variables. Namely, some SO(iV)2 states have to be projected 
out. Specifically, when N — 2n is even, the SO(2n)2 primary operators constituting the non-interacting spectrum 
(that are necessary to reconstruct the SU(2n)i Hilbert space), together with their degeneracy, are given by: 
1 x 1, {2 x }i 1 x *a++a s -,a++aj' 1 x *2A±,2A±> 1 x *2a±,2aJ ' 1 x Here the integer j runs from 1 
to n — 2, and Xj is the highest weight of the representations of SO(N), while the spinorial representations (which 
do not appear in the spectrum) have highest weight Xf . When N = 2n + 1 is odd, the spectrum is given by: 
1 x 1, {2 x $a-,a-}, 1 x *&2A s .2A a , 1 x ^ ) 2Ai.2A 1 , the integer j running from 1 to n — 1. We make the reasonable 
hypothesis that this modified SO(iV) GN model at level 2 is still fully gapped. The study of its spectrum goes 
beyond the scope of this work. 
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• AM transition: The embedding index is one in this case and one has the conformal embedding 
[51] SU(2n)i ~ Sp(2n)i x Z n , (recall that iV = 2n is even in this case) where the last piece, 
corresponding to A, stands for the parafermionic minimal model with Z n symmetry (n > 2) and 
central charge = 2(n — l)/(n + 2) which is a generalization of the Ising Z2 critical theory [52]. 
This transition was studied at length in Ref. 12, 13. In particular, it was shown that for N < 4, 
this criticality is stable, i.e. no relevant operator in the 7L n sector is allowed by symmetries that 
would spoil criticality. Once again, the quantum phase transition is associated to the spontaneous 
breaking of the Z n symmetry - with respect to the Al case, a Z2 subgroup of the center of SU(2n) 
survives that corresponds to a redefinition of the (lattice) fermion sign (ci,a — ► — Cj i(X ). 

• AM transition: The embedding index is one, and the coset model A is trivial, with = 0. 
Recall that the automorphism Q is defined by splitting the set of N fermionic chains into two 
clusters of size p and q with p + q = N. The affine symmetry g,, decomposes in three simple 
factors: su(p)i © su(?)i © u(l) (see Table I). The first two factors are generated by the affine 
currents and J^fl\, that are respectively the uniform components of the lattice spin operators 



Sa,6< P c ),a T ab c j,b and Ep< a ,6 c \,a T ab c j,b> where t(p) (respectively T^) are traceless Hermitian 



matrices with non-zero entries only in the upper- left pxp block (lower-right qxq block respectively). 
The affine current ^7°, generating the U(l) piece, is the uniform component of the relative charge 



with (j>f = (f>f L + <fif R . A very similar expression holds for the order parameter in phase ^n'- 
Ti'(r2(g)) = e V Np ^(g^') — e v N i J ^(g' 9 '). Generically, the two order parameters have 
the same symmetry and one cannot exclude coexistence. The self-dual model is described by the 
following Hamiltonian density: 



where 7~L( P ) is the SU(p) GN model at level one (similarly for 7~i( q )), so that the corresponding 
parallel degrees of freedom are fully gapped. The boson 4>f decouples and is governed by: 






(47) 



(48) 



^{^(d^ff + KfidJf)^ (49) 
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which is nothing but the Luttinger Hamiltonian, with velocity Uf = </u|- ^ 2 (pi and Luttinger 



parameter Kf = y^v F — 2irgf / yjv ¥ + 2irgf. This approach predicts a central charge c = 1 on 
the self-dual manifold. Of course, naively irrelevant operators (i.e. operators that are irrelevant 
with respect to the fully unperturbed SU(iV)i WZNW model) could be present that would spoil 
this criticality in the (j>f sector. The situation is complex in general, and we will only investigate 
here the incommensurate case, thus ignoring umklapp oscillating terms. Of course, since we start 
with a Hamiltonian local in terms of the lattice fermions, the perturbing operators can be built 
out of the fermions ^l(r)- The only SU(p) x SU(q) x U(l) candidates are of the form Ok = 
(11^=1^ ^ ) (IljLi* ^ ) +H.c, where antisymmetrization is implied in the indices a- < 
p and bj > p, and thus necessarily k < p,q. Moreover, the self-dual symmetry forces k to be an even 
integer. These operators can be expressed in terms of the SU(p)i and SU((/)i primary operators: 
Ok ~ Tr($^ Afc ) Tr($^ A *) e k 4^~^ s +H.c, where in obvious notations x k denotes the SU(p)i 
primary operator in the "p" sector. On the self-dual manifold, it is legitimate to average out the 
gapped degrees of freedom, thus yielding the following effective perturbing operator in the 4>f 
sector: 



O k ~g cosf J— — <P f ), (50) 



where g oc (Tr($fe , )) (Tr($l 9 \»)) (the average is taken in the SU(p) and SU(g) GN 
model) is a non-universal prefactor. The dominant perturbing operator is thus Oi with a scaling 
dimension (with respect to the Luttinger liquid fixed point) A2 = ^NKf/pq. We thus see that 
depending on the value of Kf the model exhibits different behaviors on the self-dual manifold: 
when Kf < pq/2N, the operator O2 is relevant, and spoils criticality: in this case, model (48) 
describes a smooth cross-over between phases ^#1 and On the other hand, if Kf > pq/2N, 

there is a true continuous phase transition, with a c = 1 criticality emerging for the degrees of 
freedom describing the relative charge on the two clusters of fermionic chains. In the latter case, it 
is instructive to consider the Hamiltonian describing the vicinity of the self-dual manifold. Writing 
the currents J±_ in terms of the SU(p)i x SU(g)i x U(l) variables yields: 

Hn = H* n + g ± Tr (g^) Tr (g^t) + H.c. (51) 

The perturbing operator is nothing but 0\. Averaging out the gapped degrees of freedom (this 
is legitimate in the limit of small g±), one concludes that the low-energy degrees of freedom close 
to the self-dual manifold, described by the bosonic field 4>f, are again governed by a sine-Gordon 
model. The scaling dimension of its cosine term Ai = NKf/pq so that a gap opens if Kf < 2pq/N 
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duality type 


criticality on the self-dual manifold 


critical degrees of freedom 


AI 

(symmetric pairing) 


no criticality (first order transition) 




AI 

(antisymmetric pairing) 


N = In < 8: Z„ parafermions 
N = 2n > 8: first order transition 


associated to the symmetry c^ a — ► e 2n / N a^ a 


AM 

(j) and q bunches) 


^ 2W : smooth crossover 
^ < Kf < 222-: c = 1 transition 
A"/ > 2£2: c=l pocket 


relative charge of the two bunches of chains 
relative charge of the two bunches of chains 



TABLE III: Different classes of quantum phases transitions between the symmetry enlarged phases, away 
from half-filling. 



as soon as g± ^ 0. If Kf > 2pq/N, the perturbation is irrelevant and one can only conclude that 
the c = 1 criticality extends on both sides of the self-dual manifold. Though this criticality is 
doomed to disappear when one departs too far from the self-dual manifold, our approach does not 
allow to put a bound to this criticality in parameter space. 

To restate our conclusions for the AM class, we find that the behavior close to the self-dual 
manifold depends on the value of the Luttinger parameter Kf, that can vary continuously. If 
Kf < pq/2N, one has a smooth crossover from phase to phase ^#r>. If pq/2N < Kf < 2pq/N, 
there is a quantum phase transition with c = 1, the self-dual model being a Luttinger liquid for 
the relative charge between the two clusters of chains. If Kf > 2pq/N, this criticality acquires a 
finite extension around the self-dual manifold. 

Table III presents a summary of the quantum phase transitions that occur between the different 
classes for general fermionic models away from half-filling. 

V. DUALITIES AT WORK: THE EXAMPLE OF THE TWO-LEG ELECTRONIC 

LADDER 

In this section we will apply the duality approach, presented in the preceding sections, to the 
study of a specific example. The novelty of this approach stems from the fact the different possible 
symmetry enlarged phases can be exhausted by means of the algebraic properties of the groups 
G and H, instead of solving numerically the one-loop RG flow and scanning the different phases. 
In addition, once a phase is identified and characterized, there is a systematic way to obtain the 
physics of the other duality-related phases. Finally, the duality approach also offers the possibility 
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to study directly phase transitions between the different symmetry enlarged phases, by investigating 
the self-dual manifolds. We now turn to a careful examination of a particular example, namely a 
generalized two- leg electronic ladder. 

A. Lattice model 

Here we apply the ideas exposed in the preceding sections to the model of two-leg Hubbard 
ladder with Hamiltonian: 

n = - Yl ( c U c i+w + H - c -) + Yl u {m°i) c \A*A,h*2 c i,hvz c iM*v ( 52 ) 

i Woo' i IjOj 

where c|» creates an electron with spin a =t,| on chain £ = 1,2, at site i. The explicit form of 
the hopping parameters t&j^'o-' and on-site couplings ^{^{{o-., } depends on the precise model under 
consideration, that can be specified by fixing the bare physical symmetry group H. 

For incommensurate filling, there is a spin-charge separation and, as already discussed in section 
II, the relevant maximal symmetry group G supported by the two-leg ladder (52) is G=SU(4) that 
mixes the four different local states {c\ ia |0) } e=l 2 -o=\ I ( nere |0) i s the vacuum state). We now 
need to specify the bare symmetry group H. Here, we have chosen to present the example of a 
model that retains three basic symmetries: first, a natural SU(2) symmetry acting in the spin 
sector; second, a Z2 symmetry that exchanges the two chains of the ladder; and finally, a U(l) 
orbital symmetry generated by N± — N2, where Ng = ^2 ia c\/> c- . is the relative charge on the two 
chains. This leads to the following symmetry breaking pattern ignoring the U(l) charge symmetry: 

G - H = SU(2) spin x U(l) orb xZ 2 x H discrete , (53) 

where H(j iscrete accounts for the remaining discrete symmetries of model (52). As discussed later, 
without affecting our conclusions regarding the number and labeling of phases, the constraint of 
the U(l) orbital symmetry can also be relaxed: in this case, and in the generic case (i.e. with a 
non- vanishing interchain hopping that is permitted in this case) , we recover in the continuum limit 
an effective U(l) symmetry group, so that the symmetry breaking pattern (53) still holds. As a 
consequence, the structure of the phase diagram, that relies on the pattern (53), will be the same. 

Given the symmetry H, we can construct the Hubbard-like Hamiltonian exhibiting this symme- 
try (for simplicity, we consider only on-site couplings) [53] . The hopping amplitudes are constrained 
to take the value twaa' = t^w^oo'- As for the terms quartic in fermions, they are obtained by 
considering the general two-fermion states at site i, c| ^cj ^ a , |0). These states decompose under 
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the group H as (1;0) © (0;0) © (0;±1), where (S;L) denotes the multiplet of spin S that carries 
orbital U(l) charge L. Denoting the 6 local two-particle states by \S, S Z ;L), they read explicitly 
as follows: 



|M;0) 


— c-^Cgj |0) 




1,-1; 0) 


= 4A io> 




|1,0;0) 


= ^ (4 T 4i 


+ c f u c 


|0;l) 


= AAi i°) 




|0;-i) 


= 4t4j 1^) 




|0;0) 


= ^ (4 T 4i 


t 



10) 



|0) , (54) 

and they allow for three different invariants under the group H, namely ^s z =o ±i 1 1> S z ' } 0) (1, S z ; 0| , 
|0, 0; 1) (0, ; 1 1 + |0, 0; — 1) (0, 0; — 1|, and |0, 0; 0) (0, 0; 0| . In this respect, we choose to parametrize 
the interacting-part of the Hamiltonian (52) as: 

Hmt = ^ [ 77 ^ ni t £ a ni t £' a ' + Jh • S;,2 + J (If) 2 , (55) 

i (for)^V') 

where rc^ = 4,&r c i,&7 is the local electronic density for species (la), S it £ = \ ^ CT(T , 4^^°%^ is 
the spin operator on chain and Tf = ^ ^ CT (n^i^ — ni,ia) is the difference of electronic densities 
on the two chains, that generates the U(l) orbital symmetry. 

The resulting Hamiltonian depends on three microscopic couplings: an (overall) Coulombic 
interaction U, a Hund coupling Jh, and an "orbital crystal field anisotropy" Jt- Playing with 
these three parameters allows to recover several limiting cases that have already been studied: 
when Jh = Jt = 0, one recovers the SU(4) Hubbard model, that has been extensively analyzed 
in recent years [42, 43, 54-56]. The case where only Hund coupling and Coulomb interaction are 
present (Jt = 0) has been studied in Refs. 17 and 57, where it has been shown that the Hund 
coupling can stabilize a (i-wave superconducting instability. When Jt = 3J///4, the symmetry is 
promoted to SU(2) S pj n xSU(2) or b, and one recovers the so-called spin-orbital model [58-62]. When 
Jt = Jh/^i the symmetry is enlarged to SO(5)xZ2, and this model has been studied in the context 
of cold-atoms physics [12, 13, 56, 63-65]. Finally, as it will be discussed later, we note that two-leg 
electronic ladders with a t± hopping process in the low-energy limit have the symmetry group (53) 
[14] so that our duality approach is also relevant to this case. 
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B. Low-energy Hamiltonian 

Let us now derive the continuous description of model (55). The procedure is standard and 
makes use of the continuum limit of the lattice fermionic operators in terms of right-and 
left-moving Dirac fermions ^£ aRtL [6, 7]: 

^ ~ * eaL (x)e-^ x + ^> i(7R (x)e +lk ^, (56) 
V«o 

with x = iao- We now introduce the SU(4)i currents which are defined from chiral fermion bilinears: 
J^r^l = ^tt^Io-b h Tfo e , a , ^e'a'R,L, T A being SU(4) generators in the fundamental representation which 
are normalized according to Tv(T A T B ) = 5 AB . The set of these 15 generators T A can be organized 
as 3 SU(2) sp j n generators T s a , a = x,y,z, that mix spin indices a, 3 SU(2) or b generators T", 
a = x,y, z, that mix orbital indices £, and 9 mixed spin-orbital generators that act simultaneously 
in spin and orbital spaces, T^, a,b = x, y, z. 

1? = ±<t*®1, T t a = ±TL®a a , T$ = ^a a ®a b , (57) 

where a a are Pauli matrices and the first (respectively second) matrix in the tensor product acts on 
spin (orbital respectively) indices. Correspondingly, we will introduce the following decomposition 
for the 15 components of the SU(4)i current: 



^SR(L) - 27r *&m(L) ( T s)la,l>a> *<Vr(l) , ^r(l) ~ 27I "^Lr(l) ( T t)ta,l>o> $ /Vh(l) 

J&M = 2 ^L( L ) { T st) w «Wh(l)- (58) 

For incommensurate filling, the low-energy Hamiltonian separates into two commuting pieces: 
TL = 7i c + TL S ([7i c ,7i s ] = 0) where the charge degrees of freedom are described by a Luttinger 
Hamiltonian: 



2 



±- (d x § c f + k c (d x e c ) 2 



(59) 



with the following Luttinger parameters for model (55): 



(QJJ J \ 

1+ao ^rJ (60) 

For generic fillings, no umklapp terms appear and the charge sector displays metallic properties in 
the Luttinger liquid universality class. All non-trivial physics corresponding to the spin degeneracy 
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is described by the spin part, i.e., TL S . Neglecting spin- velocity renormalization terms, the latter 
Hamiltonian corresponds to an anisotropic SU(4) model with marginal current-current interactions: 

5 



— 2^ {^R^R + -^L 4 ^ 4 ) + X/ 9a ^R d°AB ■Jh i 



(61) 



a=l 



where the five matrices d a encode the symmetry group (53) of the low-energy effective Hamiltonian. 
Those matrices are diagonal and they read explicitly (we choose the following ordering for the 15 
SU(4) currents: I?,I?,T%> Xt \T%): 



( 



'12 



o 6 



\ 



03 J 



J 12 



o 3 



V 







10/ 



5 



V 



(62) 



9 / 



In the continuum limit, and at first order in the lattice couplings (U,Jt, Jh)i the coupling constants 
in Eq. (61) are given by: 



91 
94 



Sir 2 
Qq 

16tt 



(2U + J t - J H ) , 92 



16tt 2 



(-4C7 + 2J t + J H ) , g 3 



8vr 2 



(2U + J t + J H ) , 



"o 



(AU - 2 J t + 3 J H ) , <?5 = (217 - 3J t ) 



(63) 



C. Duality approach 



We are now in position to apply the general duality approach to the specific model (61) to fully 
determine the nature of its spin-gapped phases. The general phase diagram of the lattice model 
(55) will then be deduced. 



1. Duality symmetries 

The Hamiltonian in the spin sector (61) takes the form of an SU(4) anisotropic current-current 
model so that we can apply the general result of section IV for the special N = 4 case. In 
this respect, we find three non-trivial duality symmetries which exhaust the possible classes for 
incommensurate filling of the general classification listed in Table I. 

A first one, to be called f2i, belongs to the AI class of symmetric pairing with gy = so (4) (see 
Table I). It is associated to the following involutive element of the center of the symmetry group 
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H: cicr — > Ri a i' a '(^,i n where the symmetric matrix R reads: R = —a y © cr y (in a good basis c a for 
the fermions, this corr 
following components: 



the fermions, this corresponds to c a — ► & a ). At the level of the SU(4)i currents, affects the 



«i = J" s f L - o, 6 = x, y, z. (64) 

This transformation is indeed a symmetry of Eq. (61) provided that the couplings are changed as 
follows: g 2 -> -g 2 and 53 -> -53. 

A second one, to be called O2; belongs to type AiT of antisymmetric pairing with gn = sp(4). 
It is associated to the following antisymmetric charge conjugation symmetry: q ct — > •^ -;' CT ' c |' ' 
where the antisymmetric matrix J (the Sp(4) metric in the C£ a basis) reads: J = io~ y © cr x . The 
SU(4)i currents affected by Vt 2 are: 

rrx rrx 

^2 : { Jl - -31 (65) 

It has the following representation on the couplings: (73 — > —g% and g$ — > —34. 

The last non-trivial duality, f2 3 , belongs to the AM class with gn = s(u(2) x u(2)) = u(l) © 
su(2) © su(2), and is associated with the orbital rotation of angle tt: q ct — > — i(— ) c^ a which is 
indeed an involutive symmetry of the bare Hamiltonian commuting with all other symmetries. It 
affects the SU(4)i currents as follows: 

J-X rrx 
th ~^ <JtL 

7 y — > — 7 y 

JtL (66) 

3 stL ~^ ~^stLi o, = x,y,z. 
This duality has the following action on the couplings: g 2 — > —g 2 and 54 — > —gn- 

Amongst those non-trivial dualities, two of them, namely Vl\ and £l 2 , are outer dualities. It is 
coherent to christen the fourth duality, namely the trivial (identity) operation, Qq. Note that the 
set of the 4 dualities {^ a }a=o,i,2,3 has the structure of the Klein four-group, or Z2 x Z2, with the 
following multiplication table for non-trivial elements: £li£l 2 = ^3, Oif^ = Q, 2 , ^l 2 Vt^ = VL\. 

Before we move on to the description of the phase diagram, it is useful to give an alternative 
form for the low-energy effective Hamiltonian. Indeed, there is a possible representation of the 
unperturbed SU(4)i model in terms of 6 free Majorana fermions, which has been widely used 
in the literature on the 2-leg ladders [16-18, 20, 22, 61, 66, 67]. This Majorana description of 
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the SU(4)i WZNW model is the quantum translation of the classical equivalence between the Lie 
algebras of SU(4) and SO(6). As we will see, this equivalence also allows for a very convenient 
representation of the possible dualities, that are realized as SO(6) dualities. 

To this end, one introduces six free massless Majorana fermions £"(r) that allow to represent 
the non-interacting part of model (61) as 

^o = -^£ - . (67) 

Under SU(2) sp i n rotations, the three Majorana £ a (a = 1,2,3) transform as a spin one, while 
the three remaining Majorana £ a (a = 4,5,6) transform as a spin one under orbital rota- 
tions of SU(2) or ij. The SU(4)i currents can then be written as fermionic bilinears. Explic- 
itly, the triplet of spin currents reads Jg h i K \ = — * 7re<lf ' c ^L(R)^L(R)' * ne triplet of orbital currents 
reads J^r R \ = — i^ eabc £ttn)£uR)> wnue the remaining 9 mixed spin-orbital SU(4)i currents read 
■^stL(R) = — ^^"(r)^^ ( a ' ^ = 1>2)3). The full interacting Hamiltonian (61) can then be conve- 
niently written as: 



H s = H s0 + 4vr 2 



l<a<6<3 l<a<3 l<a<3 

where one has introduced the energy density (or "thermal operator") in a = of the Ising 

model associated to each of the 6 free Majorana fermion theories. 

In terms of these fermionic variables, the three non-trivial dualities (64,65,66) read: 

«i = a = 4,5,6 

^2 : e L - ~e L (69) 
^ 3 :£ L a ^-£ L a a = 4, 5. 

In other words, the non-trivial dualities can be constructed from "elementary" KW dualities on the 
Ising models attached to each Majorana fermion theory. This is very peculiar to the N = 4 case, 
i.e. two-leg ladders, and this result does not generalize for N > 4. Of course, restrictions on the 
way they can be combined come from demanding that they be compatible with H-invariance. Note 
that all of these dualities, once expressed in the fermionic language, become SO(6) dualities of type 
BD7 (see Table I). By using the well-known Lie algebra isomorphisms u(l) = So(2), su(2) = so(3), 
sp(4) = so (5) and so (4) = su(2) © su(2), one can easily check that for each duality Oi,fi2,^3 
respectively, the fermionic representation (69) yields the correct invariant subspace gii, namely 
0|i = so (3) ©50 (3), 50 (5), so (2) ©so (4) respectively. 



(68) 
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2. Phase diagram 

According to the analysis of section IV, the knowledge of the different allowed dualities O a 
gives direct access to the different symmetry enlarged phases ^# a . Explicitly, the lattice order 
parameters, that will develop quasi- long-range order in phase at incommensurate filling are 
obtained from the "fundamental" CDW order parameter (Oq D ) by acting on one of the two 
involved lattice fermions with the involutive lattice symmetry associated to the duality £l a (see 
Table II). Using the dualities (69) that we have found, we can therefore readily conclude that 
model (55) with bare symmetry H (53) can support the following symmetry enlarged phases: 

• Two phases an d ^3 with coexisting SP and CDW instabilities at wave vector 2k F . Both 
CDW-SP phases are distinguished by their properties under the Z2 exchange of the two chains: ^#0 
is even while ^#3 is odd. The explicit lattice order parameters of ^#0,3 phases and their continuum 
descriptions read as follows: 

^0 : 0f DW) = O = Zia *L*&* «— e ~ l2kF3a ° Efe iia^M 

(70) 

Ufs : of^ =0 3 = EiA-Y +1 *L*toa — e-**** E to R^^A 

U ^ — e l^laK ) V eaL V iaR < e LjIv \ ) C j,ea C 3+l,£v 

• Two phases and ^#2 characterized by a superconducting pairing 0\ (O2 respectively) 
instability at wave vector k = 0, are of BCS type with s- wave symmetry (d-wave symmetry 
respectively). The two BCS phases, are distinguished by the symmetry of the pairing operators: 
0\ (BCSs) is odd under chain exchange, while O2 (BCSc?) is even. 

BCSs = Jt\ : of C) = 0i = * aL i? a6 * 6a *=- c ja R ab c j+l b = c j+1 ^c m - c j+l ll c j 2 ^ - (1 2) 

BCSd = ^f 2 : Of C) = 2 = ^>aLJab^ bR < \c ja JabC jb = C m j C m - C jiU C i2T , 

(71) 

with a = (£, a) and b = (£',&'). The BCS instability is accompanied with a generalized spin 
superfluidity, that would reveal itself for example in a rotating system put on a ring, with a 
quantization of the current in a mixed spin-orbital direction. As it has been described in section 
IV C, this result stems from the fact that Q± and O2 axe outer dualities. Moreover, if the charge 
sector is gapped by some process (for example, by considering a suitable commensurate filling away 
from half filling), one can also characterize the BCS phases by the appearance of coherence of the 
spin excitations (<5-peak in the zero temperature spin structure factor), contrarily to the CDW-SP 
phases which correspond to inner dualities. 
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The general analysis of section IV also allows for a determination of the quantum phase tran- 
sitions between the SU(4)-symmetry enlarged phases: there are as many different transitions as 
there are non-trivial dualities. Close to a quantum phase transition, the physics of the two-leg lad- 
der spin sector is captured by a minimal theory with special self-dual symmetry (which is simply 
the invariant subspace of the Z 2 grading f2 ), which is maximal in some sense: quantum phase 
transitions require in general that the bare anisotropy be large enough so that DSE cannot fully, 
but only partially, develop. What happens is sketched in the following diagram: 

/ SO(3)xSO(3) ^ 
SO(3)xU(l)xZ 2 — SO(4)xSO(2) — SU(4) , 

SO(5)xZ 2 

with the bare theory symmetry indicated on the left, and the RG flows towards low energy when 
moving to the right. In generic situations (i.e. far away from multicritical points), there is a 
hierarchy of couplings and DSE develops in several steps (in our case, two steps) as long as the 
energy cut-off is reduced. Close to a quantum phase transition this process is stopped at an 
intermediate stage, and the low-energy theory has a an intermediate symmetry that eventually 
labels the quantum phase transition. We list in Table IV the three different quantum phase 
transitions occurring in the model. 





Symmetry 


Order 


Critical degrees of freedom 


and 

^#2 *-* -^3 


SO(3)xSO(3) 


1 

2 
2 


51,54 = 95 > 
80(3)! (c = |) 5154 = 5i55 < 
so ( 6 )i ( c = 3) 5i,54 = 55 <0 


JKq «-> Jt 2 
and 

^#3 <-> ^#1 


SO(5)xZ 2 


2 
2 


Ising (c = |) 51 = 52 = 55 > 
80(6)! (c = 3) 5i= 52 = 55 <0 


.-#0 ^#3 

and 
-#2 *-* ^1 


SO(4)xSO(2) 


2 
2 


L.L. (c = 1) gi = g a > 
80(6), (c = 3) gi = 53 < 



TABLE IV: Summary of the different possible quantum phase transitions for the generalized two-leg elec- 
tronic ladder (61) for incommensurate filling. 



The next question is whether those four phases are actually realized in the lattice model (55), 
and what are the values of the microscopic couplings that will yield each of those phases. To answer 
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this question, one needs to perform, in a weak-coupling approach, a RG analysis of the continuous 
effective theory (61). The two- loop /3-function can be obtained for example in a point splitting 
regularization [24] to yield: 

91 = (1 - + 23 2 2 + S3) 

92 = 32(231 + 35) + 3334 - y ( 2 3i + 3#! + 33 + 34 + 3i) 

33 = 2#i£ 3 + 2^34 - 33(3? + 32 + 33 + 3l) ( 72 ) 

34 = 3435 + 33233 - y ( 3 32 + 3 33 + 34 + 31) 

35 = (1 - 35)(332 + 9a), 

where g = din(a/a ) , with Oq 1 the UV cutoff (which is of the order of the inverse of the lattice 
spacing) and a -1 a running RG scale, and the couplings in Eqs. (72) have been rescaled according 
to g a -> 2ng a /v F . 

One observes that the RG equations are invariant under the duality action on the couplings, 
as it should be (in fact, since dualities are exact symmetries of the continuous model (61), the all- 
order (3- function is invariant). Numerical integration of the RG equations yield the phase diagram 
presented in figure 3. The four possible symmetry enlarged phases ^#0 to ^#3 are reached within 
this model for repulsive Coulombic interaction U (see figure 3). 

DSE ensures that all bare theories sitting between the transition lines (bold lines) has a low- 
energy physics described by an adiabatic deformation of the corresponding representative (f2 a - 
twisted) SU(4) GN model. Closer to transition lines, a cross-over regime occurs (gray areas). Even 
closer to them, one finds again a universal regime described by the quantum phase transition. 
Note that the curvature of the phase transition lines visible in figure 3, is not a two- loop effect, 
but rather stems from the non-linearity of the flow equation that are already present at one-loop 
order. One notices a quite complex topography of the phase diagram, with for example an s-wave 
superconducting pocket for ferromagnetic Hund coupling. 

D. Allowing inter-leg hopping 

In this section we show how the previous discussion can also be applied to the case of a two- 
leg ladder with legs connected by transversal hopping t±, with a SU(2) spin symmetry, and a Z2 
symmetry that exchanges the two legs, but no U(l) orbital symmetry. This model has been exam- 
ined in several works, with the conclusion that there were four phases at incommensurate fillings 
[14, 18, 66, 68, 69]. We will see that this model is in fact in closely related to the SU(2) xU(l) xZ 2 
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FIG. 3: Phase diagram of model (55). The Coulombic repulsion U = O.Oli is fixed, while Hund coupling Jjj 
and orbital anisotropy J t are varied. Bold lines indicate the location of the transitions between the different 
symmetry enlarged phases whose nature is discussed in the text. The small pocket lying between points 
O and D is of type ^#3. The dashed line extending between points O and A is fully critical, with central 
charge c = 3. Thin rays correspond to fine-tuned bare theories with a larger, exact, bare symmetry H' 
with H=SU(2)xU(l)xZ 2 C ETC SU(4)=G. On those lines, the effective theory has symmetry H' at any 
stage of the RG flow, but flows at low energies to a symmetry enlarged phase with symmetry SU(4) in 
the sense of section IIIB. On point O, the bare theory is maximally symmetric, with a SU(4) invariance. 
On point A = (—8U,—2U), the continuous theory enjoys an enlarged twisted SU(4) symmetry fi2(SU(4)). 
Self-dual manifolds intersect the (Jh,Jt) plane at the discrete points B = (— If/,— |f7), C = (— 8f7, 6f7) 
and D = (0,2U), with fi x (C) = C, n 2 (B) = B and fl 3 (D) = D. 

model studied in the preceding subsections: it is simply obtained therefrom by acting with a 
duality. 

Introducing the bonding and anti-bonding modes that diagonalize the kinetic term, and expand- 
ing those modes around the corresponding two Fermi points fc F1 and k F 2, one obtains fermionic 
fields ^o-p, cr =t, I, £ = 1, 2, P = L, R. The Fermi wave vectors satisfy k F+ = k F ± + k F 2 = ^r, with 
n the electronic density, while the difference /c F _ = k F \ — k F 2 is a function of t± and n. At fillings 
and t± such that k F ^ is incommensurate, it is easy to see that for an arbitrary 2iV-fermion term 
n^Li ^I^ct'P ^e'o-'p' t° conserve lattice momentum, it has to conserve separately the combinations 
p + = Nm + N2L and p_ = A^2_r + Nil, where Np P = f{ p ^ ffP is the total number of chiral 
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fermions in mode I. Therefore, on top of conserving the SU(2) spin and the total number of elec- 
tron p + + p„ , the Hamiltonian has to conserve the combination 2j = \ (p+ — p~ ) = — ■ This 
operator is nothing but the "orbital current" , i.e. the space-like component of the Noether current 
associated with U(l) orbital symmetry. It results that the model generically enjoys in fact a larger 
symmetry SU(2)x U(l) x Z2, where the twisted orbital symmetry U(l) is generated by T t . 

Therefore, in the continuum limit one recovers exactly the model studied in the preceding 
subsections, provided one performs a duality Cl that changes the sign of the following component 
of the currents 14 : 



( 

J-X rrx 
a, ~^ Jtt, 



Jl - -Jl (73) 



'tL ^tL 
' Sth '-'stht 



3 JL ~~ ^ 3s0li CL — x,y,z. 



One thus immediately deduces that the maximal number of symmetry enlarged phases is four. 
Since (73) corresponds to an outer automorphism of su(4), one readily knows that pairing phases 
are mapped onto density-wave phases, and vice versa. Note that in terms of the six Majorana 
fermions £ a , this duality bears the very simple form: Q — ► Using the representation of Q 

on the original bonding and anti-bonding fermions: 

^Tl-^+^Il , — (74) 

one can readily identify the four phases = Cl(^ a ) of the t±- ladder, with order parameters 
O a = Cl(O a ) that read explicitly: 



1 

(7 



These phases correspond to the four phases found in generalized two-leg ladders with interchain 
hopping for incommensurate filling: coexistence of CDW and SP phases, d-wave and s-wave super- 
conducting phases (DSC and SSC respectively) and coexistence of time-reversal breaking phases 



14 Our choice is dictated by the following constraints: Cl has to preserve the SU(4) Lie algebra (it is an automorphism), 
it has to leave J$ invariant, and has to change the sign of J* . There is in fact a whole family of solutions, that 
are labelled by an angle a. Of course, our conclusions do not depend on this choice. 
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like d-density wave (DDW) and diagonal current (DC) phases (see for instance Ref. [69]). The 
explicit connection, together with lattice order parameter in terms of the original fermions c\ and 
C2 on the two legs of the ladder, is given by: 



SSC = JZ Q : A s = -O < — c i,iT c i,i| + c i,2T c i,2| 

DSC = : A d = -0 3 < — Cj^Cjpi - Cj^cj^ 

CDW+SP = Jh : O cdw = 23fie(e a ' F+x 5 2 ) <— E^R^ct^A 

5P = 4cos(A: F _f)Ke( e *+( a: +^)d 2 ) <— E^R* +1 4^ +1A + Hc " 

DDW+DC = Jg\ : O ddw = -2%m(e ik *+ x Oi) < — i E CT i„- + H - c - 

We thus see that the duality approach allows to capture the four phases of the generalized 
two-leg ladder with interchain hopping at incommensurate filling. In addition, the nature of the 
different quantum phase transitions between these phases is still determined by Table IV. Finally, 
our approach leads us to conclude that more exotic phases can only be found by (i) considering 
theories with large anisotropic bare couplings, i.e. theories that flow, in the IR limit, close to the 
quantum phase transitions described previously, or by (ii) breaking the remaining bare SU(2) spin x 
Z2 symmetry. 



VI. CONCLUSIONS 



In this paper, we have developed a general non-perturbative approach to describe spin-gapped 
phases of weakly interacting one-dimensional degenerate fermions. In the continuum limit, the low- 
energy properties of these systems are described by a WZNW CFT perturbed by marginal relevant 
current-current interaction with H invariance. At the heart of the analysis is the existence, in 
this general class of model, of emergent duality symmetries which enable one to relate different 
competing orders between themselves and to shed light on the nature of the zero-temperature 
quantum phase transitions. In particular, we have shown that these dualities can be classified and 
depend on the algebraic properties of the problem: the physical symmetry group H of the system 
and the maximal continuous symmetry group of the interaction G. For iV-component degenerate 
fermions, this maximal symmetry group is G = SU(iV) or SO(2iV) away from half-filling and 
at half-filling respectively. The duality symmetries can be identified as the involutions which 
belong to the center group C(H) Jnv of the symmetry group H of the problem so that for each 
duality, there is an involutive Z2 lattice symmetry that gives rise to it. Alternatively, based on 
the DSE phenomenon, the duality symmetries are also in one-to-one correspondence with the 
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different symmetry enlarged phases that can be supported by the generic fermionic model with 
marginal current-current interactions. Those symmetry enlarged phases are fully gapped phases, 
and all meet at the multicritical point, that is obtained by fine-tuning all interactions to zero. 
Moreover, and this should not be a surprise since we are dealing with one-dimensional models, 
those massive phases do not break spontaneously any continuous symmetries. Hence, our approach 
allows to draw a picture of a quite large class of possible generalized spin-liquid phases supported 
by coupled fermionic chains in one dimension, at least those that develop close to the multicritical 
(non-interacting) point. 

In the course of our study, there naturally appears a strong distinction between two kinds 
of gapped spin-liquid phases for incommensurate filling according to whether the duality they are 
associated to is an inner and outer automorphism. This quite mathematical distinction turns out to 
have important physical signatures. "Inner" phases, appear to display "conventional" properties, in 
the following sense: the ordering is of charge-density type, with a lattice order parameter of the form 
Cj a M a jjCj y b. On the other hand, "outer" phases display off-diagonal order, with the development 
of quasi-long-range correlations for pairing operators cj a M a ;,cj b . This superconducting instability 
is accompanied by spin superfluidity: the low-energy collective modes carry spin currents. In 
spite of the gap, the system has non- vanishing susceptibility in some spin directions, that results 
for example in a quantization of the spin current for a rotating system put on a ring. All these 
"unconventional" properties can be ultimately connected to the spontaneous breaking of a discrete 
symmetry of the bare theory: the discrete phase redefinition of the fermions (43), which it is 
tempting to connect to the usual breaking of U(l) gauge invariance in a superconducting state in 
higher dimensions. 

We can also address the issue of quantum phase transitions between the different phases, by 
means of "minimal theories" interpolating between them. The nature of the phase transition 
between the two symmetry enlarged phases and (associated to dualities £li and SI2) 

depends only on the duality 0, = so that the same set of dualities can be used to label those 

minimal models. The self-dual manifold of these minimal models exhibit in general criticality, 
that captures the nature of the quantum phase transition between the phases. We have shown 
in the specific example of two-leg electronic ladders at incommensurate filling how this approach 
allows for an immediate determination of the possible symmetry enlarged phases on the basis of 
the analysis of the bare symmetry group (of the continuous theory) and the determination of the 
different quantum phase transitions. 

As perspective, it would be interesting to study the interplay of symmetry enlargement with 
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doping - i.e. how commensurability generically affect the general picture of the phase diagram. 
In this respect, the half- filled case is very special since there is no spin-charge separation and the 
charge degrees of freedom cannot be disentangled from the spin ones due to an umklapp process. 
As it has been discussed in this paper, the relevant maximal symmetry group for N fermionic 
species at half-filling is G = SO(2iV). The detail of the analysis of the dualities for this case turns 
out to be more complicate than for the incommensurate case, i.e. with G = SU(iV), and will be 
investigated elsewhere 15 . Dualities of the kind studied in this work can have further applications 
in other contexts. For example, in degenerate quantum impurity problems where a localized spin 
is coupled to electronic spin currents, they can be used to relate different IR boundary fixed points. 
More generally, the duality symmetries can be viewed as automorphisms of the fusion algebra of 
a CFT. We focused here on WZNW models perturbed marginally by current-current interactions. 
Generalized dualities could be used to shed light on the phase diagram of other CFT's perturbed by 
relevant operators. In this respect, this case will be useful to investigate the possible classification 
of ID spin-liquid phases where the spin gap is opened by a strongly relevant perturbation. We 
hope to come back to these issues elsewhere. 
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APPENDIX A: CHARACTERIZATION OF DUALITIES 

In this Appendix, we first characterize the set of linear transformations of the currents that 
leave the Hamiltonian (5) globally invariant, by (i) showing that these transformations can be 
decomposed in the product of diagonal and chiral transformations and by (ii) showing that the chiral 
transformations, that correspond to dualities, are given by Eq. (13), i.e. they have to belong to the 
involutive part of the center of the bare symmetry group H. Then we show that for a given model, 

15 An exception occurs for N = 4, the 2-leg electronic ladder, for which the group identity SU(4)~SO(6) allows 
to factorize the SO(8) dualities into a charge part and a spin part. One can then show that with respect to 
the incommensurate case, a doubling of the number of allowed dualities occurs, yielding 8 different phases for 
the example studied in section V. It is important to realize that the triality - a remarkable, non-involutive 
automorphism of SO(8) - that was invoked in Ref.14 in the study of the half-filled 2-leg electronic ladder has 
nothing to do with the dualities studied here. Rather, it is connected to the identity SU(4)~SO(6), and helps 
understanding why N = 4 is so particular. We shall come back to this in a forthcoming publication. 
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the set of dualities coincides with the different possible inequivalent symmetry enlargements. We 
recall that Q denotes the set of allowed dualities, D is the vector space spanned by the matrices d a . 
In particular, D is the set of real (from hermiticity) symmetric (from parity) matrices that commute 
with the action of H. From renormalizability, D is closed under anticommutation: {D, D} C D. 

1. Dualities are the involutive part of the center of H 

We investigate here the general case of a transformation fi, that acts both on the currents and 
the coupling constants, such as to leave Hamiltonian (5) globally invariant. acts linearly on the 
currents, and one can write its action as uj l x u R as in section III: 

^(r) = ( w l(r)) ^ifa); ( A1 ) 

while it changes the coupling constants as g a — ► g a = £l(g a ). From the invariance of the current 
OPE's (6), it follows that both u; L and ui R must be automorphisms of G. One useful remark is that 
G-automorphisms automatically belong to 0(dim(G)), i.e. the matrices w l (r) are orthogonal. The 
covariance of the interacting part of the Hamiltonian (5) leads to an identity that generalizes Eq. 
(12): 

g a d a =g a t u R d a uo^ (A2) 

where t to R denotes the transpose of the matrix u> R . One starts by decomposing ui L x u> R = (to R x 
uj r ) o (d) x 1), with u = as the product of a diagonal and a chiral part. Considering the 

special isotropic ray g a = g%, with g^d a = 1, and multiplying (A2) by lo r on the left, and by t u> L 
on the right, one deduces that the matrix u~ l belongs to D: lu^ 1 = f a d a , with f a = Q(g^,). It 
results that is a symmetric matrix; but to being also orthogonal, one has u~ l = %j = u. We 
deduce thus 

to eB, and to 2 = 1. (A3) 

Now, we want to show that the diagonal rotation u R x u R alone leaves the Hamiltonian (5) 
globally invariant. This is equivalent to showing that t uj R d a u> R £ D for all a. Multiplying relation 
(A2) by uj r on the left and by uj r 1 on the right, one gets the following relation, valid for any 
couplings g a : 

g a oj R d a t uj R = g a d a oj. (A4) 
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Transposing this relation, and recalling that u and d a are symmetric, we deduce that uj commutes 
with d a for all a (for this, one chooses couplings gp = Q~ 1 (5 a p), so that gp = 5 a p). It results that 
d a to = \{d a ,uj} G D since {D,D} C D. We deduce then from Eq. (A4) the announced result: 

f uj R d a u R G D, (A5) 

i.e. the diagonal rotation uj r x uj r leaves the Hamiltonian (5) globally invariant. This proves our 
claim that the general transformation f2 given by Eq. (Al) can be decomposed into the product of a 
diagonal rotation, that corresponds to a global change of basis, and of a transformation oj affecting 
only one chirality sector, both of these transformations separately leaving the Hamiltonian (5) 
globally invariant. 

We now proceed to characterize the set of chiral transformations uj, that it is legitimate to term 
dualities in view of their involutive character. To study the dualities, we therefore set uj r = 1 in 
the following, with no loss of generality, so that O is a transformation of the left currents and of 
the couplings, that satisfies: 

g a d a = g a d a u. (A6) 

We already know that us belongs to D. Moreover, from Eq. (A6), fixing a and choosing gp = Q(5 a g) 
(i.e. gp = S a p), one gets d a ui = Q(5 a p)d 13 G D. Transposing this relation, and recalling that dP and 
uj are symmetric, we deduce that [u, d a ] = 0, or u>d a uj = d a , Va. It results that uj is an element of 
the physical symmetry group H. Since uj G D and [H, D] = 0, uj is in the center of H. Thus, one 
has proved that S> C C(H)| inv . 

Reciprocally, let us take some uj G C(H) . . Being involutive and orthogonal, uj has to be a 
real symmetric matrix. Furthermore, it commutes with H. This is enough to ensure that uj G D. 
To show that uj is a duality, it suffices to prove that Eq. (A6) holds. This is easily done: d a ui = 
±{uj,d a } G D. 

This ends the proof of Eq. (13). 

2. Dualities correspond to the different possible DSE's 

It is not difficult to realize that each duality fi£ f defines a different DSE. Considering the G- 
isotropic model with interacting part of the Hamiltonian H-L. = g J R J^ and acting on it with 
a duality O G £P, one gets the model with interacting part Ti^L. = g J R £l{J^) = g J R J^- 
Since Q is an automorphism of G (O G Aut(G)), the currents = 0(j7j^) are the generators of 
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the (twisted) group Gl. is thus obviously invariant under the group (Gl x GrJ^j generated 
h Y J A + J A . 

Reciprocally, let us suppose there exists some global invariance group G, isomorphic to G, that 
leaves invariant a model of the form (5), with (global) generators Q A linearly related to the original 
ones in each chirality sector, Q A = J dx (j A + J A ), with J^L\ = Then, considering the 

non- interacting part of the Hamiltonian, it follows that w L ( R ) must belong to Aut(G); considering 
the interacting part, which must be of the form g J A J A from G-invariance, leads to the conclusion 
that oj = uj r 1 u) h G D. Since ui belongs to DnAut(G), following the same line of reasoning as in A 1, 
we conclude that oj G C(H)|j = Qi . Thus, up to a global rotation uj r x uj r (affecting identically 
both chirality sectors), the generators Q A are nothing but the dual generators J dx (j A +u) AB J B ) . 
This completes the proof. 

Hence, it results that the set of dualities identifies with the different possible DSE patterns, i.e. 
they correspond to the different possibilities to glue together the two chiral invariance groups G L 
and G R in a way consistent with H-invariance. 
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